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Slow  and  steady  wins  the  race. 

Easop 


PREFACE 


The  primary  purpose  or  this  study  was  to  develop  a  computer  model  of  mutual 
coupling  in  linear  arrays  of  parallel  dipoles,  and  assess  the  impact  of  this 
coupling  on  the  side  lobe  performance  of  low  side  lobe  array  distributions  (like 
Chebyshev  amplitude  tapers).  The  secondary  purpose  of  this  study  was  to  invest¬ 
igate  simple  compensation  schemes  to  account  for  coupling.  A  fortran  computer 
program  was  written  to  model  coupling  for  any  geometry  of  parallel  dipoles,  and 
applied  to  a  variety  of  array  sizes  using  several  Chebyshev  amplitude  tapers  and 
linear  phase  tapers.  Curves  of  side  lobe  degradation  as  a  function  of  array 
size,  scan  condition  and  desired  side  lobe  level  were  produced.  Three  simple 
compensation  techniques  were  attempted. 
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ABSTRACT 


This  computer  analysis  assessed  the  impact  of  mutual  coupling  on  the  side 
lobe  performance  of  Chebyshev  linear  arrays  as  a  function  array  size,  scan  con¬ 
dition,  and  the  severity  of  Chebyshev  taper.  Array  elements  were  horizontal 
parallel  dipoles  over  a  ground  plane.  Three  simple  compensation  schemes  to 
account  for  coupling  were  investigated. 

A  fortran  computer  program,  based  upon  King-Middleton's  modified, 
zero-order,  two-term  theory  for  cylindrical  dipoles,  was  written  to  calculate 
the  generalized  impedance  matrix  for  the  linear  array,  dipole  current  distribu¬ 
tions  and  the  far-field  H-plane  pattern  in  the  presence  of  coupling.  It  is 
sufficiently  general  to  allow  the  user  to  specify,  number  of  dipoles,  dipole 
length  and  radius,  dipole  spacing,  height  above  ground  plane,  frequency,  and 
complex  generator  voltages  driving  the  array.  The  results  relate  side  lobe 
degradation  (growth)  to  array  size,  scan  condition,  and  desired  Chebyshev  side 
lobe  level. 

Additional  programs  were  written  to  reverse  the  matrix  solution,  and  solve 
for  generator  voltages,  which  after  coupling,  produce  a  base  current  distribu¬ 
tion  with  a  Chebyshev  taper.  Three  simple  compensation  schemes  then  used 
approximations  of  the  compensated  (complex)  voltages  to  drive  the  array.  The 
effectiveness  of  these  compensation  schemes  are  discussed.  The  results  of  this 
study  indicate  that  linear  arrays  greater  than  20  elements  can  probably  ignore 
mutual  coupling  for  side  lobe  considerations,  and  amplitude  only  compensation 
schemes  significantly  improve  side  lobe  performance  under  all  but  the  severest 
scan  condition  (45  deg)  and  Chebyshev  tapers  (-40  db)  considered. 


THE  EFFECT  OF  MUTUAL  COUPLING  ON  THE  SIDE  LOBE 
PERFORMANCE  OF  CHEBYSHEV  ARRAYS, 

VOu.  _ 

INTRODUCTION 

Array  analyses  that  ignore  mutual  coupling  are  based  upon  pattern  multipli¬ 
cation.  That  is,  a  common  element  pattern  assumed  to  be  applicable  to  all  array 
elements,  is  multiplied  by  an  array  factor,  which  includes  all  complex  element 
weights,  and  path  length  differences  due  to  geometry  (5,  Chap  5)(9,  Chap  4)(24, 
Chap  3).  For  arrays  of  dipoles,  this  assumption  implies  identical  current 
distributions  for  all  dipoles.  This  implication  is  not  correct.  Mutual  cou¬ 
pling  induces  parasitic  currents  on  surrounding  dipoles  due  to  to  the  excitation 
of  any  single  dipole.  The  element  pattern  of  a  driven  dipole  is  then  the  super¬ 
position  of  fields  due  to  the  driven  current,  and  all  parasitic  currents.  Since 
the  physical  and  electrical  environment  changes  across  the  array,  and  in  parti¬ 
cular  near  the  edges,  the  parasitic  currents  and  therefore  the  element  pattern 
change  as  well  (7).  The  element  pattern  cannot  be  factored  out  of  a  far-field 
expression.  This  study  will  investigate  the  impact  of  mutual  coupling  on  the 
sidelobe  performance  of  linear  arrays.  Parallel  dipole  elements  will  be  used 
due  to  their  simple  geometry  and  significant  coupling  properties.  The  utility 
of  this  can  be  extended  due  to  equivalences  between  cylindrical  dipoles,  strip 
dipole,  vertical  monopoles  over  a  ground  plane,  as  well  as  two-wire  and  wave¬ 
guide  fed  slots  (9,  Chap  7). 

The  modelling  of  mutual  coupling  has  progressed  from  analyses  that  use 
assumed  current  d; - tributions,  to  analyses  that  solve  for  actual  current  distri¬ 
butions.  The  early  work  of  Carter  (3)  used  the  induced  EMF  Method  to  calculate 
the  coupling  between  an  isolated  dipole  pair  is  an  example  of  the  former.  He 


xi 


assumed  sinusoidal  current  distributions,  and  used  the  calculated  coupling  of 
the  isolated  dipole  pair  for  any  similarly  spaced  dipoles  in  an  array.  The 
analyses  of  coupling  and  electromagnetic  fields  for  arrave  by  I4ansen  {ID 
and  King  (14,  Chap  5  and  6),  solved  simultaneously  for  all  coupling  coefficients 
between  array  elements  as  a  single  boundary  value  problem.  This  method  also 
assumed  sinusoidal  currents,  and  its  accuracy  is  limited  by  the  suitability  of 
sinusoidal  currents  to  the  dipole  geometries  under  investigation.  Later  work  by 
Chang  and  King  (4)  presented  rigorous  and  highly  accurate  solutions  for  dipole 
current  distributions  by  solving  for  coefficients  to  a  five  term  expansion  of 
dipole  current.  This  work  accurately  calculated  coupling  between  dipole 
elements  but  proved  entirely  too  cumbersome  for  application  to  array  problems. 
Later  work  by  King,  et  al  (15)  developed  a  much  simpler  two-term  theory  for 
current  distibutions,  and  applied  this  theory  to  linear  arrays  of  thin,  parallel 
Hi pol ps.  Thi<=  simpler  theory  achieved  excellent  agreement  with  Mack's  precise 
measurements  (18)  of  dipole  impedances  and  current  distributions. 

Section  1  outlines  fundamental  electromagnetic  equations  for  the  cylindrical 
dipole  antenna.  Section  2  and  Appendices  Al  and  A2  present  King-Middleton' s 
Modified,  Zero-Order,  P/ro-Term  Theory.  This  theory  is  used,  in  Section  3,  to 
model  the  coupling  of  parallel,  horizontal  dipoles  over  a  ground  plane.  A  vari¬ 
ety  of  array  sizes  (10,  20,  30,  and  40  elements)  employing  numerous  Chebyshev 
amplitude  (15,  20,  30,  and  40  dB)  and  linear  phase  (0,  15,  30,  and  45  degrees  of 
scan)  tapers.  Section  4  discusses  the  impact  of  mutual  coupling  on  the  sidelobe 
performance  of  the  Chebyshev  tapers,  as  a  function  of  array  size  and  scan 
condition.  Section  5  employs  three  simple  compensation  methods  to  correct  for 
the  coupling,  and  discusses  their  effectiveness.  Appendices  B  and  C  contain 
model  results,  in  the  form  of  far-field  plots  and  active,  current  •.■•eights,  for 
uncompensated  and  compensated  arrays  respectively.  Volume  II  of  this  report 
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1.0  FUNDAMENTAL  ELECTROMAGNETIC  THEORY 

The  fundamental  characteristics  of  cylindrical  antennas,  the  electromagne¬ 
tic  fields  ana  drivinq  point  impeaancss,  are  completely  determined  by  the 
distribution  of  currents  on  the  antenna.  The  first  section  of  this  report  will 
review  the  pertinent  electromagnetic  equations  for  both  isolated  and  coupled 
cylindrical  antennas  (5) (9) (15) (24) .  Impedance  and  far-field  expressions  will 
be  derived,  assuming  an  idealized  "slice"  generator. 

Fields  problems  are  usually  expressed  by  linear  integral  equations  derived 
from  Maxwell's  equations  and  specific  boundary  conditions.  Maxwell's  equations 
can  be  stated  as  follows  for  homogeneous  mediums. 


V  x  E  =  -joc^H  (l-2a) 

V  •  E  =  p/e  ( l-2b) 

V  X  H  =  jooeE  +  J  ( l-2c) 

7*5=0  (l-2d) 


J  is  assumed  to  be  conduction  current,  related  to  the  electric  field  by  the 
conductivity,  a  .  That  is, 

J  =  <jE  (1-3) 

For  good  conductors,  J  is  constrained  to  be  surface  currents,  as  interior  elec¬ 
tric  fields  are  zero.  The  volumetric  charge  density,  p  ,  is  also  constrained  to 
be  surface  charge  for  good  conductors.  The  permeability  and  permitivity  are  // 
and  e  ,  respectively.  All  field  properties  (  E,  H,  J  and  p  )  are  infact  both 
spatial  (x,  v  and  z)  and  temporal  (t)  functions.  The  time  dependence,  e;,wt:  ,  is 
assumed  for  ail  field  properties  and  suppressed  throughout  this  report. 
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1 . 1  HELMHOLTZ  EQUATIONS 

Since  E  and  H  are  coupled  in  Maxwell's  equations,  the  equations  must  be 
solved  simultaneously  to  determine  rbe  fields  from  *  given  current  distribution, 
J  .  Therefore,  we  introduce  scalar  and  vector  potential  functions,  <t>  and  A  ,  to 
decouple  and  simplify  the  solution  for  the  fields  t24:§  1.3).  To  satisfy 
( 1— 2d ) ,  we  define  the  magnetic  field,  H  ,  as  a  curl  of  the  magnetic  vector 
potential  function,  A  .  That  is, 


H 


1 

U 


r  ' 

7  X  A 


(l-4a) 


By  integrating  (l-2a),  we  can  form 

E  =  - jccA  -  V<fr  ( 1— 4b) 

where  the  electric  scalar  potential  function,  4>  ,  is  yet  to  be  fully  determined. 
To  fully  determine  $  ,  we  impose  the  Lorentz  condition  upon  A  and  <)>  ,  that  is 

V  •  A  =  (l-4c) 

substituting  (l-4a  and  c)  into  (l-2c)  we  form  the  Helmholtz  equation  for  vector 
potential  (or  the  vector  wave  equation), 

V2  A  +  B2A  =  —pj  ( l-5a ) 


where 

|32  =  (p  pt 

Likewise,  substituting  (l-4b)  and  (l-4c)  into  (l-2b),  we  form  the  Helmholtz 
equation  for  scalar  potential  (or  scalar  wave  equation), 

V2  +  62  <j>  =  -p/e  .  ( l-5b) 
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For  good  conductors,  the  distributions,  J  and  o  ,  become  the  surface 
distributions,  J  and  o3  .  The  particular  integrals  for  (l-5a  and  b)  can  now  be 
expressed  in  terms  ot  the  given  surface  distributions  and  the  free  space  Green- s 
function  for  outwardly  radiating  spherical  waves.  That  is, 


A(r) 


(  r 


J3  (r'  ) 

/ 


s 


e 


-  j  PR 


4IIR 


ds' 


and 


4>(  r )  =  - 
c 


p.  (r' ) 


3  0R 


4 IIR 


ds' 


( l-6a ) 


( l-6b) 


where 


R  -  |r  -  r'  | 

which  is  the  distance  between  the  field  point  and  the  source  point.  Lastly,  the 
law  of  conservation  of  charge  is  imposed  upon  J  and  p  .  Namely, 

7  •  j  =  —  j  cop  .  (l-6c) 

At  the  interface  between  two  media,  regions  1  and  2,  the  following  boundary 
conditions  must  be  satisfied: 


-  E2  )  X  ft  =  0 

( l-7a) 

ft  X 

(El  -  E2  )  =  ps/e 

( l-7b) 

ft  X 

<sx  -  s, )  -  J 

( l-7c ) 

ft  X 

( Hx  -  IL  )  =0 

( l-7d) 
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where  n  is  the  unit  normal  vector  to  the  surface,  pointing  from  region  1  to 
region  2. 

1.2  THIN-WIRE  APPROXIMATIONS 

Thin-wire  approximations  are  not  limited  in  application  to  thin-wire  struc¬ 
tures,  as  general  computer  codes  have  proven  themselves  useful  for  modelling 
arbitrary  conducting  surfaces  with  thin-wire  grids  (20:§  2.7).  If  the  wire 
radius,  a  ,  is  small  relative  to  wavelength  and  wire  junction  separations  are 
long  relative  to  radius,  then  current  flow  is  axial  and  uniformly  distributed 
around  the  wire  circumference.  Further,  all  currents  and  charges  can  then  be 
approximated  by  equivalent  line  sources  located  along  the  wire  axis.  In  other 
words, 


1(1)  =  2IIaJg  (1) 
q(l)  -  2nap8(l)  . 

From  (l-6c),  we  can  restate  the  law  of  conservation  of  charge. 


( l-8a) 
{ l-8b) 


j  3  Id) 

q(  1 )  =  -  - 

(o  91 


( l-8c ) 


Finally,  the  thin-wire  solutions  to  the  Helmholtz  equations  can  be  written  in 
terms  of  line  integrals, 


A(  r)  =  p 


Kl' ) 


-  j@R 


4  riR 


dl' 


wire 


(l-8d) 


<P(  r) 


J 

we 


9  1* 


1(1' ) 


31 ' 


j  0  R 


dr 


4IIR 


wire 


(l-8e) 
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These  integrals  express  the  potential  functions  m  terms  of  a  known  source  {cur¬ 
rent  distributions  in  this  case).  The  boundary  condition  to  be  satisfied  for 
ihin-wires,  from  (1-t?).  4  <= 


1  •  ( E,  -  E,  )  =  0 


(l-8f) 


where  1  is  a  unit  vector  along  the  wire  axis.  For  perfect  conductors,  E.  s  0  , 
and  E,  can  be  expressed  as  a  sum  of  incident  and  scattered  fields  such  that 


E,  (1' )  +  Es  (1' )  =  0 


( l-8g) 


1.3  ISOLATED  DIPOLE 

The  simple,  realizable  transmitting  system  in  Figure  1  will  be  replaced  by 
the  idealized  linear  antenna  of  Figure  2.  In  this  analysis,  it  is  useful  and 
convenient  to  replace  the  actual  transmitter  and  transmission  line  with  a 
"slice"  generator  (described  in  equations  I-9a  and  b). 


T  ransmitter 
RF  Source 


Transmission  Line 


z  —  +  h 


z  =  0 


FIGURE  1.  REALIZABLE  LINEAR  TFAMSMI TIER 
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ing  wire  (or  tube),  of  length  2h  and  radius  a,  with  hemisherical  end  caps.  This 
antenna  is  bisected  by  a  "slice"  voltage  generator  which  maintains  a  uniform 
voltage  step  between  the  arms  of  the  antenna.  Both  the  density  of  current,  J, 
and  the  volume  density  of  charge,  p,  are  exactly  cero  within  a  cerfect 
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conductor.  In  general,  the  boundary  conditions  for  normal  and  tangential  elec¬ 
tric  (E)  and  magnetic  (H)  fields,  at  the  interface  between  two  regions 
(14:12-13) .  can  be  expressed 


A1 

X  Ex 

+ 

A2 

X 

E: 

=  0 

(l-10a) 

A1 

X«i 

+ 

A2 

X 

*2 

=  -K1H2) 

(l-10b) 

A1 

•  A 

+ 

A2 

• 

*2 

■  -(psl+Ps2)/£o 

( l-10c) 

A1 

+ 

A2 

• 

a. 

4 

=  0 

( 1-1 Od) 

where  h1  is  the  outward  unit  vector  normal  to  the  region  1,  Ja  i  is  the  density 
of  surface  current,  and  psi  is  the  density  of  surface  charge  at  the  interface. 
If  region  1  is  defined  to  be  the  antenna  (perfect  conductor),  and  region  2  is 
defined  to  be  air,  the  boundary  conditions  reduce  to, 


A2 

X  e2 

=  0 

(1-lla) 

A2 

x^ 

- 

(1-llb) 

A2 

•  *2 

=  -p  ,  /e 

si'  0 

(1-llc) 

A2 

'  % 

=  0 

(1-lld) 

Equation  (1-lla)  states  that  the  tangential  E-field  at  the  surface  of  the 
antenna  is  continuous  and  zero.  Equation  (1-llb)  states  that  the  tangential 
H-field  is  discontinuous  and  proportional  to  the  surface  density  of  current  on 
the  antenna.  Likewise,  (1-llc)  states  that  the  normal  E-field  is  discontinuous 
and  proportional  to  the  density  of  surface  charge  on  the  antenna.  Equation 
(1-lld)  states  that  the  normal  H-field  is  continuous  and  zero  at  the  surface. 

For  thin,  Z-oriented,  perfectly  conducting,  cylindrical  antennas,  where  all  vec¬ 
tor  potential  is  also  Z-oriented,  the  Laplacians  of  ( l-5a  and  b)  reduce  to, 
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( l-12a) 


(  3‘/  Sz*  +  0‘  )A„  =  0 
(  32/3z2  +  0o2  )*  =  0  .  ( l-12b) 

The  particular  integrals  of  (l-12a  anu  D;  can  oe  derived  for  Z-onented  antennas 
from  (l-8d  and  e)  and  they  are 


A 


Z 


4ji 


h 

/ 

lz( 2'  ) 

✓ 

-h 


-je0R 

e  dz' 

R 


1 


4ite 


h 

q(z'  ) 


where  the  total  axial  current  is  from  (l-8a), 


Iz  (z)  =  2naJz (z) 
and 


( l-13a) 


(1-1 3b) 


(l-13c) 


0  -  2it/X 

o  o 

R2  -  (z-z')2  +  a2 

q(z)  is  the  distribution  of  charge 

From  (l-4a  to  c),  we  can  totally  determine  the  fields  (E  and  H)  from  the 
vector  potential,  Az .  That  is,  for  Z-oriented  dipoles  (in  cylindrical  coordin¬ 
ates  ) , 

1  9A 

Hp  = - - 

P  9<t> 

-1  3A 

H.  -  —  7^ 

iU,P  9p 
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H  =>  0 

Z 


-jw  f  3A„ 


P  |32  Up3z 

O 


-ju  (1  3  A 


Ip  3<J>3z 


■j«  r9  A,  2 

-  -  +  BA 

00  ^  ^  02 


For  rotationally  symmetric  cylindrical  antennas,  all  partial  differentials  with 
respect  to  vanish,  therefore 


H  -  H  -  0 

P  z 


(l-14a) 


-1  3A 

H(  .--1 
P0  3p 


( 1— 14b) 


E+  -  0 


(l-14c) 


-j«  A 


p  ft  bp3z 


(l-14d) 


-jw  (d  A 

E  -  -  - -  +  B2A 

ft  l  3z2  02 


(l-14e) 


Translating  (l-14a  through  e)  into  spherical  coordinates,  we  have 


H,  -  *.  -  0 
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-1  3A 


( l-15a) 


H* 

/«  3p 


E*  “  0 


Ef  -  Ez  COS©  +  Ep  SIN© 


COS© 


a2  a,  p2Az 

- -  +  e2A  +  SIN©  — - 

2  o  z 

3z  I3p3z 


( l-15b) 


Ez  SIN©  +  Ep  COSO 


COSO  -  -  SIN© 

Up3zJ 


—  +  0  2A 

2  =  * 


( 1-1 5c) 


At  large  distances  from  the  antenna  (r2  ))  h2  &  (|3or)2  ))  1)  t  the  field 
equations  reduce  to  the  simpler  radiation  or  far-field  expressions  (15:6-7), 
where  uniform  wave  front  approximations  result  in  only  transverse  (but  ortho¬ 
gonal)  E  and  H  fields.  That  is 


y  -  -  *e;/c 


(l-16a) 


y  -  ©e; 


3«^0  f  -je0R 

-  SIN©  I  (z' )e  dz' 

4n  z  R- 


(l-16b) 


From  Figure  3,  the  distance  R,  can  be  expressed  via  the  cosine  law 


R2  =  r2  +  (z')2  -  2rz'C0S© 


(1-17) 
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In  the  far-field,  where  r2  ))  (z')2  ,  the  linear  terras  (in  z')  of  a  binomial 
expansion  of  (1-17)  are  sufficiently  accurate.  Therefore, 

R  a  r  -  z'  COS©  .  (1-18) 

This  approximation  has  the  useful  property  of  separating  the  slowly  varying 
amplitude  dependence  (in  r),  from  the  phase  dependence  in  (l-16b).  That  is, 

F  (9,8  h)  (l-19a) 

O  0 


where 

h 

6  SIN©  f  j(3  z '  COS© 

F  (0,6  h)  =  — -  I(z')e  °  dz' 

0  0  _  z 

2 

-h 


( l-19b) 
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The  "field  factor",  Fo'0,f3oh)  ,  embodies  all  directional  characteristics  of  the 
electric  field  of  a  Z-oriented  dipole  of  length  2h.  Likewise,  the  magnetic 
field,  in  the  radiation  zone,  is  at  right-  e.r.giec  b’ih  rr^portionaJ  to  the  elec¬ 
tric  field.  From  ( 1— 16a )  we  see  that  the  H-field  shares  the  same  directional 
characteristics,  or 


J 

B ♦  3  7 

2rt 


(  -je0r- 


F  (0,6  h) 

o  o 


:i-i9c) 


To  formulate  expressions  for  driving  point  impedances,  we  must  use  the 
concept  of  the  complex  Poynting  vector,  and  the  EMF  method  (15:9-11).  The  com¬ 
plex  Poynting  vector  is  defined  as  follows, 


S  s 


E  X  B* 


o 


EXH 

2 


(l-20a) 


It  is  convenient  to  think  of  S  as  a  vector  field  indicating  the  direction  and 
magnitude  of  power  flow.  Now  the  EMF  method  solves  for  time  average  power 
transmitted  by  an  antenna  (more  properly  by  a  transmitting  system),  by  taking 
the  real  part  of  the  integral  of  the  normal  component  of  S,  over  a  surface  that 
completely  encloses  a  transmitting  system  (usually  a  great  sphere). 


where 


(l-20b) 


ft  is  an  outwardly  pointing  unit  vector  of  the  surface 

Zt  is  a  surface  that  completely  encloses  the  transmitting  system. 
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If  the  maximum  current  value,  Im  (  @  2=0  for  h  <  a/4  or  @  z  =  n-X/4  for  h  >  X/4, 
given  sinusoidal  currents),  is  related  to  the  total  power  radiated  through  the 
surface  cf  a  qreat  spnere,  the  concept  of  radiation  resistance,  R' ,  can  oe 
defined 

P  s  1  I  Im  |2Re  .  (1-21) 

2 

The  value  of  R*  is  not,  in  general,  the  same  as  the  driving  point  resistance, 

Ro .  Only  under  the  conditions  of  a  sufficiently  thin  (  a/X  <  10' 5  ),  resonant 
dipole  ((3oh-»rt/2)  does  the  radiation  resistance  approximate  some  of  the  driving 
point  circuit  properties  on  the  antenna  (15:14-16). 


It  is  more  common,  when  computing  the  circuit  properties  of  linear 
antennas,  to  integrate  S  over  the  cylindrical  surface  of  the  antenna  (15:22-23). 
That  is  for  a  Z-oriented  dipole. 


where 


h 


-rca 


E  H*  dz 
z  * 


-h 


(1-22) 


Sp  is  the  p  directed  component  of  S 

Ez  is  the  E-field  on  the  surface  of  the  conductor 

is  the  conjugate  of  the  H-field  on  the  surface  of  the  conductor. 


Setting  (1-22)  equal  to  common  circuit  equations  for  the  "slice"  voltage 
generator  terminals. 


13 


h 


1/2 


=  1/2  V  I*  =  -rca 

0  O 


E 

z  * 


-n 


dz 


Solving  for  the  driving  point  impedance,  z  ,  we  have 


Z 


h 


E  H. 

Z  * 

-h 


dz 


(1-2 3a) 


( l-23b) 


The  H-field  on  the  cylindrical  surface  of  the  conductor  cam  be  expressed  by 
plugging  (l-13c,  l-14a  and  b)  into  1-llb  and  solving  for  H* , 


therefore 


(1-230 


Z  = 

O 


E  I*dz 

Z  Z 


-h 


This  approach  to  the  driving  point  impedance  has  ignored  the  boundary  con¬ 
dition  stated  in  (l-8f)  (15:25-29),  which  states  the  tangential  E-field  vanishes 
at  the  surface  of  a  perfect  conductor.  Since  (l-8f)  must  be  satisfied,  some 
impressed  E-field,  E* ,  supported  by  some  continuous  distribution  of  voltage 
generators  is  also  implied.  The  impressed  E-field  exactly  cancels  the  E-field 
in  ( 1— 14e ) ,  such  that 

E  +  E9  3  0  :.  E9  3  -E 

z  z  z  z 
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Such  a  continuous  distribution  of  generators  is  not  conducive  to  the  derivation 
of  driving  point  properties,  as  there  is  no  unique  driving  point.  The  center 
driven  antenna  with  vanishing  tanaem-i  (except  in  the  gap)  is  there¬ 

fore  fundamentally  different  from  an  antenna  with  tangential  E-fields  expressed 
by  (l-14e  and  l-13a)  due  to  assumed  sinusoidal  currents.  The  currents  and 
fields  are  not  generally  the  same,  however  under  certain  circumstances,  some 
properties  are  comparable.  For  resonant  dipoles,  the  current  amplitude  and 
phase  distributions,  of  the  sinusoidal  current  model,  have  proven  (via  measure¬ 
ment)  to  be  reasonably  accurate.  Characteristics  that  depend  wholely  upon  the 
current  distribution  (I  ,  Z  ,  &  H. )  must  also  be  reasonably  accurate.  Since  the 
E-field  in  the  radiation  zone  is  proportional  to  the  magnetic  field,  it  is  com¬ 
parable  as  well.  This  cannot  be  said  for  the  full  wavelength  (anti resonant) 
dipole. 

For  near  resonant  dipoles,  even  though  currents,  near  zone  3-fields  and  all 
far  zone  fields  are  reasonably  approximated  by  sinusoidal  currents,  measured 
charge  distributions  are  very  disimilar  [15:27],  at  the  center  and  ends  of  the 
antenna,  from  the  cosinusoidal  charge  distribution  required  by  the  condition  of 
conservation  of  charge.  Since  the  radial  E-field,  Ep ,  is  proportional  to  the 
charge  per  unit  length  along  the  antenna,  near  zone  E-fields  will  not  be  compar¬ 
able. 


1.4  COUPLED  DIPOLES 

If  it  were  not  for  the  coupling  between  antennas,  the  driving  point  proper¬ 
ties  of  coupled  antennas  would  be  identical  to  the  isolated  antenna,  and  the 
far-fields  would  be  a  straight  forward  superposition  of  the  fields  from  each 
antenna  taken  separately.  Unfortunately,  the  antennas  do  couple,  so  the  cur- 
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rents  cn  one  antenna  are  dependent  cn  the  excitation  ot  all  antennas  (to  some 
degree),  and  therefore  the  fields  cannot  be  formulated  from  the  superposition  of 
isolated  antennas. 

To  compute  the  seif  and  mutual  impedances  of  the  coupled  dipoles,  the  EMF 
method  is  again  employed.  This  time  however,  the  Z-oriented  E-field  at  the  sur¬ 
face  of  one  antenna,  is  the  sum  of  fields  from  that  antenna's  own  currents,  and 
the  currents  of  all  other  antennas.  Consider  the  case  of  two  coupled  dipoles, 
in  Figure  4. 


FIGURE  4.  COUPLED  DIPOLES 


where 


( l-24b) 


PK  is  the  total,  time  averaged  power  crossing  though  the  surface  of 

the  Kth  dipole 

Sp  is  the  radially  directed  Poynting  vector 

Ezk.  is  the  E-field  at  the  surface  of  the  Kth  dipole  due  to  currents  on 

the  ith  dipole 

H+ki  is  the  H-field  at  the  surface  of  the  Kth  dipole  due  to  currents  on 
the  ith  dipole. 

Mote  that  the  only  H-field  encircling  the  Kth  dipole  is  due  to  currents  on  the 
Kth  dipole  (15:29)  (see  l-23c),  (l-24a)  simplifies  to 

Pk  ’  j  J  i  (E,I“I“) dz  '  (1~24c) 

-h 


Equating  (l-24c)  to  the  power  expression  for  the  driving  point, 


1  1 1  .  1 2 Z  .  =  11  Cl.  i*Z.. 'I 

1  °k  °k  I  i  (  1  k  lkJ 


where 


(l-25a) 


( l-25b) 
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and 


Zok  is  the  drivinq  ooint  impedance  of  the  Kth  dipole 

Zkk  is  the  seif  impedance  of  cue  KLii  dipoie 

Zki  is  the  mutual  impedance  between  the  ith  and  Kth  dipoles 

Iok  is  the  complex,  driving  point  current  amplitude  of  the  Kth  dipole 

Ik  is  the  complex  current  amplitude  of  the  Kth  dipole  due  to  the 

excitation  on  the  Kth  dipole 

Just  as  in  the  isolated  dipole  case,  the  sum  of  all  E-fields  (E  )  must 
vanish  on  the  surface  of  all  dipoles,  except  in  the  gap  region.  Note  again, 
however,  that  the  charge  distributions,  and  the  radial  components  of  the  near 
zone  E-fields,  are  not  well  approximated  even  for  near  resonant  dipoles.  This 
has  somewhat  greater  impact  on  coupled  antennas  (15:44-45),  because  near  zone 
fields  become  important  in  determining  parasitic  (coupled)  currents.  The  far- 
field  is  a  superposition  of  far-field  from  all  currents  (base  driven  and  parasi¬ 
tic).  King-Middle ton's  two  and  three-term  theories  address  this  problem  in 
particular,  in  that  they  solve  for  the  currents  on  each  dipole  (allowing  those 
currents  to  vary  across  the  coupled  elements)  and  then  derive  the  driving  point 
and  far-field  properties  from  these  more  accurate  currents. 

Even  though  the  assumed  sinusoidal  currents  enable  the  approximation  of 
driving  point  self  and  mutual  impedances  (from  1-27),  they  are  more  central  to 
the  computation  of  far-fields  in  conventional  array  theory.  Apply  (l-16b)  to 
the  coupled  dipole  case  in  Figure  4.  The  far-field  is  now  a  sum  of  integrals 
(one  over  each  element's  current  distribution). 


18 


(( 

El  =  - -  SIN©  I. 

4rt  . 


where  (see  Figure  5] 


-jMi  !" 

z '  )e  dz '  -  i 

R. 


-^,R:  .  ,  ! 
z' )e  az' 


(1-26) 


Rj  =  r-z'COS0 


Rj  =  r-z'COS0  -  bC0S4>J 

R1  =  R2  =  r  in  the  denominators 


in  the  exponential  term 


Dipole  1 


/^R2  ~  r-bcos  <j> 
*  x 


Dipole  2 


boos  b 


FIGURE  5.  FAR-FIELD  APPROXIMATION  FOR  R1  &  R2 


Let  each  element  current  include  an  arbitrary,  complex  amplitude  (current 
weight),  such  that 


!.z(z')  «  Iio  e  ‘  Iiz  (z' 


Then  (1-26)  could  be  rewritten, 
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(1-27) 


By  assuming  all  current  distributions,  l.z(z')  ,  are  identical,  we  can  separate 
Eg  into  a  field  factor,  Fm(©,0oh),  which  embodies  all  directional  properties 
attributable  to  the  current  distribution  (assumed  to  be  applicable  to  all  ele¬ 
ments),  and  an  array  factor,  A{0,<j>),  which  accounts  for  all  time  delays  (or 
equivalently  phasing)  due  to  geometric  placement  of  each  element  (eJpobcos+)  and 
the  complex  current  weights  (I^e3*1)  (see  equation  l-28a). 
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m  o 


( l-28a ) 
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6,  SIN© 
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( l-28b) 


For  a  linear  array  of  N  elements,  the  array  factor  could  be  extended. 


A(©,$) 


N 
=  Z 
i=l 


/ 

I. 

1  O 


j(si+e0 


( i-l)bC0S<J>)' 
/ 


(1-29) 


In  actuality,  the  current  distributions  on  coupled  antennas  vary  with  the 
geometric  and  electromagnetic  environment  in  an  array  [14:72-73].  Once  again, 
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Kina-Middletcn's  two  and  three-term  expansions  for  oyiindricai  antennas  will 
allow  for  more  accurate  far-field  ccmputaticns.  This  accuracy  is  gained  at  the 
expense  of  evaluating  an  ’Pr?gr ai  ■ n  • L~i~ 1  tor  everv  ccuoiec  antenna. 

1.5  RUMSEY  REACTION  INTEGRAL 

The  Rumsey  Reaction  Technique  is  an  alternative  physical  interpretation, 
that  is  useful  in  this  (and  many  other)  electromagnetics  analyses.  Reaction  is 
basically  a  measure  of  coupling  between  one  source  and  another  (20:§  2.13)(24:S 
7.4).  Consider  a  wire  structure  (Figure  6)  enclosed  by  a  surface  S  ,  in  the 
impressed  fields,  E;  and  H.  ,  due  to  remote  sources,  J.  and  M.  . 


Total  Raids 

E  =  Ei  i-  Es 
u  -  u;  > 


FIGURE  6.  SCATTERER  IN  IMPRESSED  FIELDS 


The  total  fields  are  the  sum  of  incident  fields,  and  scattered  fields  due  to  the 
wire  scatter.  From  Schelkunoff rs  surface  equivalence  theorem,  the  wire  scatter 
can  be  replaced  by  equivalent  surface  sources,  Js  and  Ms  (see  Figure  7).  The 
equivalent  sources  are  defined  by  the  curl  expressions. 


J  =  ft  x  H  ( l-30a) 

S 

M  =  E  x  ft  .  ( l-30b) 
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E  =  Ei 
H  =  Hi 


[ero  -ields 


*  Eauivelant 

Js.  Ms  Surface  Sources 


Remote 

Sources 


FIGURE  7 .  EQUIVALENT  SURFACE  SOURCES 

Now  add  to  this  problem  test  sources,  Jt  and  ,  that  are  located  within  S  (see 

Figure  8).  The  total  fields  now  include  a  term  for  the  fields  due  to  the  test 
sources. 


Jt.Mt 

Ef.Ht 


Total  Fields 

E  =  Ei  +  £3  +  Et 
H  =  Hi  +  Hs  r  Ht 


.  Equivelant 

Js.  Ms  Surface  Sources 


Remote 
Sources  t 


FIGURE  8.  ADDITION  OF  TEST  SOURCES 


The  reaction  between  the  external  sources,  J  J  H  and  M_  ,  and  the  test 
sources,  J.  and  M.  ,  can  be  expressed  in  the  manner  of  Rumsey, 
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<i+S,t>  -  E. +E  -J  -  H  +H  • M  dr  =  0 

i  •»  t  i  -  r 


:  l-31a] 


By  reciprocity,  we  can  also  write 


<i+s,t>  s  <t,i+s>  = 


E.  •  J.  +J 

t  IS 


'  _  1  ' 

-H,  •  M.  +M  dT 

t  13 

v  /  . 


l-31b) 


For  the  purposes  of  this  development,  assume  only  electric  test  sources,  only 
impressed  electric  sources  and  a  perfectly  conducting  scatterer.  Therefore, 


Mj.  =■  0 


M.  »  0 


M  -  0 

S 


Equation  (l-31b)  reduces  to, 


Et  •  J.  +  Js  dx  =  0 


(l-32a) 


Expanding  (l-32a)  by  integrating  over  electric  sources  separately,  we  have 


Et  -  J.  dx  +  Et  -Js  ds  =  0  . 


( 1— 32b ) 


Applying  reciprocity  to  the  integral  over  the  impressed  sources, 


E.  -  J  dt  +  E/j  ds  =  0 

1C  C  3 


:i-32c) 
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For  a  perfectly  conducting  scatterers,  the  test  sources,  J„  ,  will  be 
constrained  to  the  surface  also,  therefore, 


f 

(  ' 

f 

r 

r 

E  •  J 

ds  =  - 

E  • J 

i  t 

t  3 

> 

V  V 

> 

> 

k  > 

t  s 


( l-32d) 


Equation  (l-32d)  is  the  Rumsey  integral  equation,  which  can  be  solved  for  Ja  in 
terms  of  known  E  ,  and  chosen  J  -*  E.  .  For  thin-wire  structures  equation 
(l-32d)  can  be  reduced  to  line  integrals.  The  Rumsey  Reaction  Concept  is  dis¬ 
cussed  here  because  it  is  central  to  Richmond's  Thin-wire  Method  of  Moments 
Code.  That  established  code  will  be  used  later  as  a  check  this  anthor's  imple¬ 
mentation  of  King-Middleton' s  Two-Term  Theory  in  this  study. 
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:.0  KING-MIDDLETOM'S  .".EDIFIED,  ZERO-ORDER,  CEO-IEEM  THEORY 

King  (14)  ,  and  King,  Hack,  and  Sandler  13),  after  long  and  careful 
~orT"p?r'scn  ~z  linear  an  re  nr.  a  analyses  uicr.  ...easu.  whence,  r^oucsca  an  ^rararr”e 
solution  to  Hallen's  integral  equation  for  currents  on  a  cylindrical  dipole. 

This  solution  assumed  a  three-term  form  for  the  dipole  current  distribution,  and 
then  solved  a  boundary  value  problem  to  determine  the  coefficients  for  each 
term.  King,  et  al,  demonstrated  for  near  resonant  isolated  dipoles,  and  for 
linear  arrays  of  parallel,  near  resonant  dipoles,  that  a  condensed  two-term  for¬ 
mulation,  requiring  a  lumped  capacitive  correction  to  self  admittances,  is 
sufficiently  accurate.  This  section  presents  the  relevant  two-term  equations 
(derived  in  Appendix  Al),  for  isolated  and  coupled  dipoles.  King-Middleton' s 
Modified,  Two-Term  Theory  (as  implemented  by  this  author)  will  be  compared  with 
J.H.  Richmond's  well  known  Thin-Wire  Method  of  Moments  Code,  developed  at  Ohio 
State  Universities  Electroscience  Laboratory  (21) (22). 

2.1  THE  ISOLATED  DIPOLE 

Consider  again  the  idealized  dipole  in  Figure  9.  The  Hallen  Integral  Equa¬ 
tion  relates  the  magnetic  vector  potential  due  to  currents  along  a  dipole,  to  a 
general  solution  to  the  Helmholtz  equation  for  magnetic  vector  potential. 
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'  -j4n 

f  Vo 

I, ( z' )K( z, z' )dz'  =  - 

c  cose  z  +  —  sins  \z\ 

'  2  3  J 

-n 


where  K(z,z')  is  the  free  space  Green's  function.  In  equation  (2-1),  C  has  yet 
to  be  determined.  King-Middleton  actually  used  a  difference  integral  equation, 
that  enabled  them  to  solve  for  all  coefficients. 


h 


/ 

f  ' 

— j  4  rt 

f  v  i 

o 

I#(*') 

K(z,z' )-K(h,z' ) 

dz'  =■  - 

U  F  ( z )  —  nz(z) 

J 

n  cose  h 

l  2  J 

where  Vo  is  the  complex  gap  voltage,  no  is  the  intrinsic  impedance  of  free 
space,  and 


h 

-jh  ( 

U  -  - 1  I  (z' )K(h,z' )dz'  ( 2-2b) 

4n  J 

-h 

F  (z)  =  COS0  z  -  COS0  h  ( 2-2c) 

oz  o  o 

Mob(z)  -  SIN60(h-|z|)  .  (2-2d) 

To  reduce  (2-2a)  to  an  algebraic  expression,  King-Middleton  expanded  the 
integral  into  four  integrals,  over  the  real  and  imaginary  parts  of  the  Kernal, 
K(*,z').  Their  study  of  the  behavior  of  each  integrand  enabled  them  to  approxi¬ 
mate  each  integral  by  an  algebraic  term  (involving  an  expansion  coefficient). 
This  analysis  led  King-Middleton  to  a  three-term  current  form  for  isolated 
antennas,  which  collapsed  to  a  two-term  form  for  near  resonant  dipoles.  By 
inserting  this  two-term  form  into  (2-2a),  and  grouping  coefficients  of  like 
terms,  King-Middleton  could  solve  for  the  coefficients  of  the  two  current  terms. 
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They  had  thereby  determined  the  current  distribution  along  the  dipole.  Their 
two-term  theory  resulted  in, 


I 

Z 


(z) 


3  2tl  Vo 

%  ^dR 


Go  (z)  +  T'  F  ( z ) 

O  Z  O! 


(2-3a) 


where 

(*v  -  SINeoh>  SEC  60h  +  *du  SIN0o h  -  jtdI 
T'  -  -  (2-3b) 

\p  -  'PJ„  cose  h 

u  d  U  o 

Goz(z)  =  SIN0jz|  -  SIN0oh  ( 2-3c) 

and  <j/dR  ,  \j/v  ,  ,  ^d(J  ,  and  ^dI  are  defined  in  Appendix  Al  (equations  A-25a 

to  h).  Given  the  current  distribution  in  (2-3a)(in  terms  of  the  dipole  geometry 
and  the  gap  voltage)  one  can  solve  directly  for  the  driving  point  properties  (at 
z  -  0)  ,  or  integrate  over  the  current  to  compute  the  farfields.  The  expression 
in  (2-3a)  includes  some  algebraic  and  trigonometric  manipulation  to  remain 
determinate  at  0oh  =  n/2  (see  Appendix  Al). 

In  Figure  10,  the  author 'a  two-term  calculations  are  compared  with 
King-Middleton's,  two-term  and  three-term  theories,  as  well  as  Mack's  precise 
measurements,  as  a  function  of  dipole  length  (in  wavelengths)  and  a  fixed  dipole 
radius  ( a/X)  .  The  conductances  are  seen  to  agree  very  well,  while  the  sub- 
ceptances  are  seen  to  be  somewhat  in  error.  Figure  11  compares  the  author's 
two-term  calculations  with  the  output  of  Richmond's  thin-wire  code  for  the  same 
dipole  geometry.  Richmond's  Thin-Wire  Code  is  derived  from  Rumsey  Reaction 
Concepts  applied  to  thin-wire  structures.  The  test  sources  and  remote  sources 
are  chosen  to  be  overlapping  sinusoidal  dipoles  and  a  generalized  impedance 
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matrix  is  produced  for  all  sources  (21).  This  is  equivalent  to  a  Method  of 
Moments  procedure  using  sinusoidal  basis  functions  and  Galerkin's  method 
•  21 )  ( 12 ) .  Again,  the  ^-thcr  'c  tvo-Lerr  -l culat’ compare  favorably. 
King-Middleton  modified  their  two-term  theory  with  a  lumped  capacitive  (subcep- 
tive)  correction  at  the  gap.  King  (15)  supports  this  modification  by  discussing 
the  inability  of  an  ideal  "slice"  generator  (infinitesimal  discontinuity  in 
scalar  potential  at  the  dipole  center)  to  account  for  terminal  zone  effects  of 
an  actual  dipole  and  driving  circuitry.  King  has  successfully  matched  measured 
'apparent'  admittances  by  a  lumped  subceptive  correction.  Their  correction  was 
chosen  to  match  measurements  for  0  h  =  ri  .  This  author  chooses  to  correct  his 

O 

two-term  implementation  at  0  h  =  1.44  ,  for  resonant  dipoles  which  shall  be 
used  in  later  array  analysis.  Figure  12  shows  that  a  lumped  correction  of  -2.27 
mMhos  should  be  used. 

By  plotting  the  current  in  (2-3a)  against  Mack's  measurements  for  a  half 
wavelength  dipole,  we  note  good  agreement  (see  Figure  13).  We  expect  the 
far-fields  from  these  similar  currents  to  agree  equally  well. 
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TH  KING-MI DDLETCN  THEORY  AND  MACK'S 


TWO-TERM  THEORY 


Dipole  Current  Distribution 
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THE  COUPLE  DIPOLE 
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Consider  the  coupled  dipole  pair  in  Figure  14.  .or  coupled  dipoles,  che 
— o<?“J  r'-;cs::cial  at  the  surface  of  one  ciccle  -.3  cue  •'  ■-  ■rr ~i!--u  .1  a_. 

dipoles. 


P(x...r: 


i  e 
hs. 
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Diooie  2 


FIGURE  14.  COUPLED  DIPOLE  PAIR 


Therefore  (2-2a)  would  be  written  for  dipole  1, 
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f  4 

f  \  ‘ 

Iiz(z') 

Ki: (z,z' )  -  KL1 (h,z'  M  +  I,z  (z' ) 
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/ 

-h 


-j  4n  f  V 

-  u,fo  (z)  +  —  n  (z) 

n  cose  h  l  2 

O  0 


( 2-4a ) 


and  for  dipole  2, 
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-h 
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K22(z,z') 

\ 

-  K.  (h,z'  ' 
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-j  4rt  '  V 

-  U.  F  ( z )  +  — 1  M  (  z  ) 

^  O  Z  ^  0  2 

n  COS  6  h  ,  2 


1 2-4b) 


where  Ux  ,  U2  ,  K11(z,z')  =  K22(z,z')  and  K,  2  ( z ,  z ' )  3  K21(z,z')  are  defined 
in  Appendix  Al.  Assuming  a  two-term  current  form, 


I.  (z)  =  jA.  M  (z)  +  B.  F  (z) 

12  J  1  C  Z  10Z 


(2-5) 


we  can  substitute  into  (2-4a  and  b)  and  apply  the  appropriate  integral  reduction 
rules  from  Appendix  Al.  Grouping  coefficients  of  like  terms  we  can  solve  for 
A.  's  and  Bi 's.  In  matrix  form  we  have, 


n  *  cose  h  V  , 

o  a  R  o  o  2 


(2-6a) 


( 2-6b) 


vrtiere  the  matrices  [  ]  and  (  <|>v  ]  are  defined  in  Appendix  Al.  The  matrix 

relationship  between  dipole  base  currents,  Io  ,  and  gap  voltages,  V  ,  can  be 
expressed  in  terms  of  a  generalized  admittance  matrix,  (  Y  ]  .  That  is, 


Ie  -  t  Y  I  V 


(2-7a) 


where 


K  F  (O)  +  I  M  ( 0 ) 

J  0  Z  O  Z 


( 2-7b) 


n0  ^dR  coseoh 


and  (  I  ]  is  the  identity  matrix. 


34 


Equation  ( 2— 7b )  is  indeterminate  at  0oh  =  n/2  (half  wavelength  dipoles). 
To  avoid  this  difficulty  a  special  form  can  be  formulated,  producing 


TU 
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F (o) 

o  z 


G  (o) 

o  z 


2  it 


*dR 


( 2-8a ) 


where 


cose  h 


(2-8b) 


Tv 


SIN6o  h 


(2-8c) 


For  this  special  form  we  can  rewrite  the  two-term  current  for  the  ith  dipole, 
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where 


( 2-9a ) 


A.' 
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2r 

*dR 


V  . 

O  1 


B'  = 


( 2— 9b) 


( 2-9c) 


This  matrix  formulation  is  equally  valid  for  any  linear  array  of  parallel 
dipoles.  The  diagonal  elements  of  (  Y  ]  are  self  admittances  and  would  be 
modified  with  the  capacitive  correction  of  Section  2.1.  The  off-diagonal  ele¬ 
ments  are  mutual  admittances  between  dipoles.  The  inverse  of  [  Y  ]  would  be 
the  impedance  matrix,  [  z  ]  . 
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To  impedance  match  the  ith  voltage  generator  with  the  ith  dipole,  its  internal 
impedance  should  be  chosen  as  the  conjugate  of  the  ith  dipoles  self  impedance. 

Tn  actuality,  the  active  impedance  (driving  point  impe-'-nce  g.  'en  drwi-.g 
voltages  across  the  array,  commonly  uniform;  should  be  matched.  As  active 
impedances  vary  with  frequency  and  scan  condition  (7),  that  match  is  normally 
made  for  center  frequency  and  broadside  conditions  (endfire  condition  for 
endfire  arrays). 

Figure  15  plots  self  and  mutual  admittances  for  King-Middleton  modified, 
two-term  theory,  Richmond's  thin-wire  code,  Mack's  measurements,  and  the 
author's  calculations  for  a  pair  of  full  wavelength  dipoles,  as  a  function  of 
the  dipole  separation,  b/X  .  The  self  admittances,  as  in  the  isolated  dipole 
case,  include  a  capacitive  correction  to  the  self  subceptance. 

Table  I  compares  the  active  (driven)  admittances  and  impedances  across  a  10 
element  linear  array  of  full  wavelength  dipoles,  as  computed  by  King-Middleton 
Two-Term  Theory  (unmodified)  (16),  OSU  Thin-Wire  Code  (less  King-Middleton 
Subceptive  Correction)  and  the  author's  Two-Term  calculations  (unmodified). 

Table  entries  are  in  good  agreement.  The  bar  chart  in  Figure  16  compares  the 
magnitude  of  the  dipole  active  base  currents  (normalized  to  the  center  elements) 
for  a  uniformly  illuminated  10  element  array.  Again  all  three  methods  agree 
very  well. 

To  compute  the  far-fields  of  coupled  dipoles  one  must  superpose  the  contri¬ 
bution  of  each  current  term  from  each  dipole.  For  the  dipole  pair  (N-2)  ,  we 
can  write 
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(2-10) 


where  G  (0, 3  h)  and  F  (0,(3  h)  are  as  defined  in  Appendix  Al.  In  this  work, 

mo  mo 

the  coordinate  system  is  referenced  to  the  gap  of  dipole  1,  so  is  the 

separation  between  the  gap  of  the  ith  dipole  and  the  gap  of  dipole  1. 
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Theory  Comparison  for  Dipole  Pair 


FIGURE  15.  COMPARISON  OF  COUPLED  DIPOLE  ADMI' 


10  ELEMENT  BROADSIDE  ARRAY 
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TABLE  I.  COMPARISON  OF  ACTIVE  IMPEDANCES  OF  UNIFORMLY  DRIVEN,  10  ELEMENT  LINEAR  ARRAY 


DRIVING  POINT  CURRENTS 


FIGURE  16.  COMPARISON  OF  ACTIVE  BASE  CURRENT  FOR 

UNIFORMLY  DRIVEN,  10  ELEMENT  LINEAR  ARRAY 


3.0  LINEAR  ARRAY  OF  PARALLEL  DIPOLES 

To  extend  the  theory  of  section  two  to  the  analysis  of  larger  linear  array 
geometries,  ore  starts  vrith  a  rralti -p^rt  tatri r  model  of  »n  isolated  array  of 
tuned  dipoles,  then  models  a  ground  plane  via  image  dipoles.  The  generalized 
admittance  matrix  is  computed  using  King-Middleton's,  modified  two-term  theory, 
for  all  dipoles  real  and  imaginary.  The  array  is  element-by-element  impedance 
matched,  for  the  uniformly  illuminated  broadside  condition,  by  choosing  the  vol¬ 
tage  generator's  internal  impedance  based  upon  each  elements  active  impedance 
under  broadside,  uniformly  driven  conditions.  Lastly  various  Dolph-Chebyshev 
amplitude  tapers  (for  sidelobe  control)  and  linear  phase  tapers  (for  scanning 
control)  will  be  applied  via  the  voltage  generators.  Array  size,  in  conjunction 
with  amplitude  and  phase  tapers,  will  be  varied  to  investigate  coupling  effects 
on  low-side  lobe  array  design. 


3.1  LINEAR  ARRAY  MODEL 

Figure  17  illustrates  a  geometry  for  a  linear  array  of  N  parallel  dipoles. 
Figure  18  illustrates  the  2N  port  matrix  model  of  this  array.  For  this  study, 
the  array  geometry  will  be 


1.44  -»  h/X  =  0.2291831 

(3-la) 

a/X  =  0.007022 

(3-lb) 

d/X  =  0.5 

( 3-lc) 
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Dipole  N 


FIGURE  IT.  ISOLATED  LINEAR  ARRAY  ON  M  PARALLEL  DIPOLES 


The  generalized  admittance  matrix,  [  Y  ]  ,  would  be  computed  using 
King-Middleton  Theory.  The  addition  of  a  ground  plane,  1/4  wavelength  behind 
t-hp  array  (  s/A  =  0.5  )  ,  is  accomp).  i  shed  using  image  dipoles  --dth  i:*vc-rt2'1 
generator  voltages  ( 9 : §  2.5)  (see  Figures  19  and  20). 

The  voltages,  Vq.  ,  in  Figure  20  imply  ideal  voltage  generators  with  zero 
internal  impedances.  A  constant  power,  voltage  source,  which  includes  an  inter¬ 
nal  impedance,  is  more  realistic.  To  include  constant  power  sources  in  the 
array  model,  inspect  Figure  21  for  the  equivalent  circuit  of  the  ith  dipole  and 
its  voltage  source.  A  common  voltage  divider  relationship  exists  between  the 
generator  voltage,  Voi  ,  and  the  gap  voltage,  VV  .  This  relationship  can  be 
expressed, 


V'. 

o  i 
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o  i 


2. 


i  i 


=  V 


gi 


gi 


gi 


(3-2) 


The  gap  voltages  from  (3-2)  can  now  be  applied  to  the  admittance  matrix  as  in 
Figure  22. 
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N  +  2 


r 


FIGURE  19 .  LINEAR  ARRAY  OF  HORIZONTAL  DIPOLES  OVER  A  GROUND 
PLANE  MODELLED  BY  IMAGE  ARRAY  (AT  Y  -  -S) 


Dipole  12  N  N+1  N  +  2  2N 


FIGURE  20.  4N  PORT  MATRIX  MODEL  OF  LINEAR  ARRAY  OVER  GROUND  PLANE 


44 


FIGURE  21.  EQUIVALENT  CIRCUIT  FOR  iTH  DIPOLE  AND 
CONSTANT  POWER,  VOLTAGE  SOURCE 
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FIGURE  22.  4N  PORT  MATRIX  MODEL  OF  LINEAR  ARRAY  IN  FIGURE  19, 

APPLYING  NEW  GAP  VOLTAGES  TO  IMPEDANCE  TO  MATRIX  [  Y  ] 

3.2  IMPEDANCE  MATCHING 

Due  to  mutual  coupling,  currents  are  induced  cn  a  dipole  by  ail  voltage 
sources  in  the  array.  Impedance  matching  the  ith  --oltage  generator  to  the  self 
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impedance  of  the  ith  dipole  is  therefore  insufficient.  Impedance  matching  to 
the  "active"  impedance  of  the  ith  dipoie  is  better.  Active  impedances  are  com¬ 
monly  measured  (or  calculated)  for  uniform  voltage  distributions  •■•ith  no  phase 
slope  (no  scanning).  The  active  impedances  can  be  calculated  from  (A-44b)  as 
follows, 


where 


,  active  base  currents 


(  Y  ]  =*  generalized  admittance  matrix  from  (A-4b) 


( 3-3a) 


( 3-3b) 


( 3-3c) 


V  -  r  1  ]  =  uniform  voltage  distribution  for  real  ( 3— 3d ) 

.  and  image  arrays 

1 

-i* 

.  -i  . 


and 


Y  .  =  i  . /V  =  active  admittance  (3-4a) 

a  i  a  i  /  u  i 

Zai  =  1/Ya  i  =  active  impedance  .  ( 3— 4b) 

Choosing  the  ith  voltage  generator's  internal  impedance,  Z_i  ,  to  be  the 
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conjugate  of  the  ith  dipole's  active  impedance,  maximizes  the  power  transferred 
to  (radiated  by)  the  dipole,  and  therefore  impedance  matches  the  voltage  genera¬ 
tor  for  uniformly  illuminated,  broadside  conditions.  Assuming  tr.e  array  vi11 
scan  some  volume  about  broadside,  this  is  also  a  reasonable  match  over  expected 
scan  conditions.  Note  that  both  the  currents  and  driving  voltages  for  image 
dipoles  are  negative,  therefore  the  active  impedance  of  these  dipoles  and  their 
images  are  equal.  Substituting  Y*.  for  Yg.  in  (3-2)  incorporates  an 
element-by-element  active  impedance  match  into  the  4N  port  matrix  model,  illu¬ 
strated  in  Figure  22. 

3.3  CHEBYSHEV  AMPLITUDE  TAPER 

The  generator  voltages,  Voi  ,  in  Figure  18,  are  complex  array  weights  that 
can  be  varied  to  control  features  of  the  array's  farfield  pattern. 
Dolph-Chebyshev  amplitude  tapers  are  real  (amplitude  only)  tapers  developed  by 
Dolph  (8)  and  pursued  by  Brown  and  Sharpe  (2)  as  a  realizable  manner  to  optimize 
either  beamwidth  given  a  sidelobe  constraint,  or  sidelobe  levels  given  a 
beamwidth  constraint.  The  far-fields  of  Chebyshev  line  sources  (and  linear 
arrays)  are  characterized  by  uniform  sidelobe  levels  and  the  narrowest  possible 
beamwidth  given  that.  sidelobe  level  (see  Figure  23). 

There  are  approximate  methods  (10) (1)  for  synthesizing  Chebyshev  array 
weights,  however  the  methods  result  in  approximate  Chebyshev  patterns,  having 
non-uniform  sidelobe  structure.  To  obtain  true  Chebyshev  patterns,  cumbersome 
perturbation  techniques  are  applied  to  alter  the  zeroes  of  the  approximate  pat¬ 
tern.  Sharpe  and  Brown  (2)  have  already  accomplished  this  procedure  for  a 
variety  of  array  sizes  and  sidelobe  levels.  In  this  analysis,  their  tabulated 
array  weights  will  be  used  as  voltage  generator  weights  (see  Figure  24  to  27). 
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deal  Chebvshev  Tapers 


FIGURE  23.  IDEAL  CHEBYSHEV  PA' 


FIGURE  25.  IDEAL  CHEBYSHEV  WEIGHTS,  20 


Chebyshev  Array  Weights 


ELEMENT  # 


3.4  LINEAR  PHASE  TAPER 

The  Chebyshev  pattern,  synthesized  in  Section  3.3,  can  be  scanned  in  $  by 
»<Wry?  =>  linear  phase  slope  (taper)  across  the  array.  Consider  ‘•he  ;;y  n 

Figure  28,  in  which  the  dipole  elements  are  oriented  perpendicular  to  the  page. 
If  the  phasing  of  elements,  5,  ,  is  advanced  to  the  physical  length  (in  wave¬ 
lengths)  between  the  plane  wavefront  (A-B)  and  the  ith  element,  the  array 
pattern  would  be  scanned  in  the  direction.  For  the  ith  element  the  phase 
advance  required  would  be, 

S1 

COS<f>  °  -  (3-5a) 

(i-l)d 

or 


6.  »  (i-l)d  CCS<t>  .  (3-5b) 

To  apply  such  a  phase  taper  to  a  Chebyshev  amplitude  taper,  we  advance  each 
element  weight  according  to  (3-5b),  that  is 


O  i 


(  V 

\ 

0  1 

CHEBYSHEV; 

f  v 

'i 

0  1 

\ 

CHEBYSHEV, 

6, 


COS  ((  i-1  )d  COSfj  +  jSIN  j'U-Dd  COS<j>j 


3-6) 


The  weights  in  (3-6)  will  scan  the  Chebyshev  pattern  to  the  angle  £  •  In 
reality,  the  mutual  coupling  and  the  pattern  will  vary  (distort)  with  scan.  It 
is  the  challenge  of  the  antenna  designer  to  choose  array  geometries  and  tapers 
suitably,  or  to  limit  scan  volume  such  that  the  resultant  array  patterns  are 
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FIGURE  28.  LINEAR  PHASE  SLOPE 


still  useful.  In  the  analysis,  the  main  beamwidth  will  be  seen  to  increase  and 
become  asymmetric,  but  remain  useful.  See  Figure  29  for  example  scanning  of  40 
db  Chebyshev  linear  array  of  40  isotropic  sources.  Of  more  interest  in  this 
work  is  the  growth  in  sidelobe  levels  as  scan  angle  increases. 
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Array  Pattern 


FIGURE  29.  IDEAL  SCANNING  OF  LINEAR  ARRAY  OF 

ISOTROPIC  ELEMENTS  (NO  GROUND  PLANE) 


4.0  MODEL  RESULTS 

After  running  four  linear  array  sizes  (10,  20,  30,  and  40  elements)  through 
t-he  coupling  model  applying  four  Chebyshev  amplitude  tap-"  ir,  20,  30,  und  40 
db)  at  four  scan  conditions  (0,  15,  30,  and  45  degrees)  for  each  array,  two 
expected  trends  are  evident.  First,  the  smaller  the  array  the  greater  the 
impact  of  mutual  coupling  on  the  far-field  patterns.  Second,  the  greater  the 
scan  (steeper  the  linear  phase  taper)  the  more  impact  on  sidelobe  level.  This 
is  due  in  some  portion  to  the  ground  planes  impact  on  mainlobe  field  strength  as 
it  is  scanned.  However,  that  portion  due  to  the  ground  plane  is  limited  to  ~  1 
db  over  a  45  degree  scan  angle  (see  Figures  30  and  31).  These  trends  are  illu¬ 
strated  in  Figures  32  through  35.  The  X-axis  is  the  desired  (designed) 

Chebyshev  taper  applied  to  the  array  as  real  voltages.  The  Y-axis  is  the  resul¬ 
tant,  maximum  sidelobe  level  minus  the  desired  sidelobe  level.  The  Y-axis  is 
therefore  a  measure  of  sidelobe  growth  due  to  coupling.  Each  figure  has  four 
curves  by  array  size.  Figure  32  through  35  are  for  0  to  45  degree  scan  condi¬ 
tions  respectively. 

The  following  general  comments  are  based  upon  a  fairly  arbitrary  3  db 
growth  standard.  The  sidelobe  growth  for  30  and  40  element  arrays,  due  to 
mutual  coupling,  is  minimal.  The  20  element  array  has  significant  sidelobe 
growth  only  under  the  45  degree  scan  condition,  and  then,  only  for  the  40  db 
taper.  The  10  element  array  has  significant  growth  under  all  but  broadside  scan 
conditions,  but  normally  only  for  40  db  Chebyshev  tapers  (both  the  40  db  and  30 
db  tapers  at  the  45  degree  scan  condition  exceed  the  3  db  standard). 
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Ideal  Chebvshev  Tapers 


FIGURE  31.  IDEAL  CHEBYSHEV  TAPERS  (15  DB  TAPER 
10  ELEMENTS,  45  DEG  WITH/WITHOUT  GRl 


SIDELOBE  DEGRADATION 


DESIRED  CHEBYSHEV  TAPER 


SIDELOBE  DEGRADATION 


DESIRED  CHEBYSHEV  TAPER 


SIDELOBE  DEGRADATION 
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Figures  36  through  39  are  bar  charts  of  ideal  Chebyshev  voltage  weights 
versus  the  magnitude  of  resultant  current  weights  due  to  coupling.  The  impact 
of  coupling  manifests  itself  primarily  in  the  5  edge  elements  for  ■-eve'-e  taoei  s 
(30  and  40  db,  see  Figure  39).  For  less  severe  tapers  (15  and  20  db,  see  Figure 
38),  eight  edge  elements  are  significantly  impacted  by  coupling.  The  less 
severe  Chebyshev  tapers  are  characterized  by  radically  large  edge  element 
weights,  and  are  therefore  more  impacted  by  coupling. 

Figures  40  through  47  demonstrate  that  nearly  all  far-field  patterns  are 
well  behaved  even  accounting  for  coupling.  The  only  notable  exception  is  Figure 
43.  The  40  db  Chebyshev  taper  for  small  arrays  translates  to  a  relatively  large 
main  beamwidth.  At  large  scanning  angles  beam  is  further  broadened  and 
distorted  and  vestiges  of  this  beam  have  been  scanned  out  of  visible  space. 

This  condition  is  contributing  to  the  increased  sidelobe  growth  evident  in 
Figure  43.  The  significantly  larger  growth  for  40  db  tapers  and  small  arrays 
(Figures  30  through  33,  10  element  curve)  far  exceeds  the  growth  evident  for 
even  the  45  degree  scan  condition  for  an  idealized  array  of  isotropic  radiators 
over  a  ground  (see  Figure  30). 

Figures  44  through  47  demonstrate  that  large  array  patterns  are  well 
behaved  for  all  tapers  and  scan  conditions  investigated.  This  would  indicate 
that  mutual  coupling  maybe  ignored,  for  sidelobe  considerations,  for  the  30  and 
40  element  linear  arrays. 
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FIGURE  36.  CHEBYSHEV  ARRAY  WEIGHTS  (10  ELEMENT 
0  DEGREE  SCAN,  15  AND  20  DB  TAPER) 
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FIGURE  37.  CHEBYSHEV  ARRAY  WEIGHTS  (10  ELEMEN 
0  DEGREE  SCAN,  30  AND  40  DB  TAPER 
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15  AND  20  DB  TAPER) 


FIGURE  39.  CHEBYSHEV  ARRAY  WEIGHTS  (40  ELEMENT 
0  DEGREE  SCAN,  30  AND  40  DB  TAPER) 


Chebyshev  Patterns 


FIGURE  40.  SCANNED  CHEBYSHEV  PATTERNS  (15  DB  TAPER 
10  ELEMENTS,  0  AND  15  DEGREE  SCAN) 


Scanned  Chebyshev  Patterns 


FIGURE  41.  SCALED  CHEBYSHEV  PATTERNS  (15  DB  TAPER 
10  ELEMENTS,  30  AND  45  DEGREE  SCAN) 


Chebyshev  Patterns 


SHEV  PATTERNS  (40  DB  TAPER, 
0  AND  15  DEGREE  SCAN) 


Scanned  Chebyshev  Patterns 


FIGURE  43.  SCANNED  CHEBYSHEV  PATTERNS  (40  DB  TAPER 
10  ELEMENTS,  30  AND  45  DEGREE  SCAN) 


Scanned  Chebyshev  Patterns 
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Scanned  Chebyshev  Patterns 


FIGURE  45.  SCANNED  CHEBYSHEV  PATTERNS  (15  DB  TAPER 
40  ELEMENTS,  30  AND  45  DEGREE  SCAN) 
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CHEBYSHEV  PATTERNS  (40  DB  TAPER 
24TS,  0  AND  15  DEGREE  SCAN) 


Scanned  Chebyshev  Patterns 


FIGURE  47.  SCANNED  CHEBYSHEV  PATTERNS  (40  DB  TAPER 
40  ELEMENTS,  30  AND  45  DEGREE  SCAN) 


5.0  COMPENSATION  METHODS  AND  RESULTS 

The  impact  of  mutual  coupling,  evident  in  Figures  32-35,  may  preclude  the 
dec-fan  of  1  ow  sidf=  lobe  arrays  for  given  array  size,  scan  conditio-*  and 
width  specifications  (6)  (17) (23).  The  Chebyshev  voltage  distributions  driving 
the  linear  arrays  in  Section  4.0,  do  not  translate  into  Chebyshev  base  current 
distributions  due  to  mutual  coupling.  This  section  will  discuss  three  simple 
compensation  schemes  to  correct  for  mutual  coupling,  and  then  apply  these 
schemes  to  the  linear  arrays  under  investigation. 

5 . 1  COMPENSATION  SCHEMES 

Ml  three  compensation  schemes  revolve  around  equations  3-3a  through  d. 

That  is,  having  modelled  coupling  in  the  formulation  of  the  generalized  admit¬ 
tance  matrix,  [  Y  ]  (see  equation  2-8a),  one  can  simply  reverse  the  matrix  pro¬ 
blem  (13) (23)  to  solve  for  the  complex  generator  voltages,  VV,  required  to 
achieve  a  known  Chebyshev  base  current  distribution,  Ioi  .  One  must  solve  the 
matrix  equation 

V'  -  [  Y  ]-1  I  (5-1) 

The  resultant  generator  voltages,  V'  ,  are  as  previously  stated,  complex. 
Such  an  excitation  may  not  be  realizable.  To  alleviate  this  concern,  three  sim¬ 
ple  compensation  schemes,  that  approximate  the  excitation,  ,  will  be  used  in 
this  study.  The  first  scheme,  'amplitude  only',  uses  only  the  magnitude  of  the 
complex  voltages,  V'  .  The  second  scheme,  '5  deg  phase',  adds  some  phase  con¬ 
trol  by  using  both  magnitude  and  phase  of  V'.  ,  however  phase  is  'discretized' 
to  the  nearest  multiple  of  5  degrees.  The  third  scheme,  '2.5  deg  phase',  merely 
increases  the  resolution  of  the  phase  increment  to  2.5  degrees.  In  all  three 
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,  are  used  as 


schemes,  the  approximate  compensated  generator  voltages,  7' 
driving  voltages  for  the  arrays  of  interest. 


5 . 2  COMPENSATION  RESULTS 

The  need  for  compensation  schemes  is  inversely  proportional  to  array  size, 
and  proportional  to  scan  angle  and  Chebyshev  taper  (see  Figures  32-35).  The 
effectiveness  of  the  simple  compensation  schemes  presented  here  is  just  the 
opposite  (see  Figures  48,  53-55,  60,  63  and  66).  These  simple  schemes  are 
particularly  effective  at  broadside,  regardless  of  array  size  and  taper  (see 
Figures  48-52,  60-62,  67  and  68).  This  trend  agrees  with  Herd's  experimental 
work  with  even  smaller  linear  arrays  of  dipoles  (13). 

Based  on  the  cases  modelled  in  this  study,  linear  arrays  of  20  dipoles  and 
larger  do  not  benefit  enough  from  these  compensation  schemes  (Figures  60,  63  and 
66)  to  justify  the  design  effort.  Large  scan  angles  for  the  20  element  array 
(Figure  63)  are  an  exception.  The  smaller  10  element  array  benefited  consider¬ 
ably  from  these  schemes,  however  (see  Figures  48-60)  amplitude  only  compensation 
is  sufficient.  In  general,  the  5  degree  phase  compensation  proved  considerably 
less  effective  than  the  amplitude  only  scheme.  This  is  dramatically  true  for 
severe  Chebyshev  tapers.  For  the  severe  tapers  on  small  arrays,  the  2.5  degree 
phase  compensation  was  required  to  achieve  the  desired  side  lobe  levels.  In 
general,  the  2.5  degree  phase  compensation  did  no  better  (and  often  a  little 
worse)  than  amplitude  only.  Small  arrays  with  large  scan  angles  and  severe 
tapers  were  not  helped  significantly  by  any  of  these  compensation  schemes  (see 
Figures  58  and  59).  This  is  more  probably  due  to  scanning  some  vestige  of  a 
large  mainbeam  out  of  visible  space. 
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COMPENSATED  SIDE  LOBES 


DESIRED  CHP3YSHEV  TAPER 


Compensation  Schemes 


FIGURE  49.  COMPARE  UNCOMPENSATED  AND  AMPLITUDE  ONLY  COMPENSATION 
0  DEG  SCAN,  15  DB  TAPER,  10  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  50.  COMPARE  5  DEG  &  2.5  DEG  COMPENSATION 
0  DEG  SCAN,  15  DB  TAPER,  10  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  51.  COMPARE  UNCOMPENSATED  AND  AMPLITUDE  QMLY  COMPENSATION; 
0  DEG  SCAN,  40  DB  TAPER,  10  ELEMENTS 
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FIGURE  56.  COMPARE  UNCOMPENSATED  AND  AMPLITUDE  ONLY  COMPENSATION; 
45  DEG  SCAN,  15  DB  TAPER,  10  ELEMENTS 


FIGURE  57.  COMPARE  5  DEG  &  2.5  DEG  COMPENSATION  S< 
45  DEG  SCAN,  15  DB  TAPER,  10  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  58.  COMPARE  UNCOMPENSATED  AND  AMPLITUDE  ONLY  COMPENSATION 
45  DEG  SCAN,  40  DB  TAPER,  10  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  59.  COMPARE  5  DEG  &  2.5  DEG  COMPENSATION  Si 
45  DEG  SCAN,  40  DB  TAPER,  10  ELEMENTS 
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FIGURE  61.  COMPARE  UNCOMPENSATED  AND  AMPLITUDE  ONLY  COMPENSATION 
0  DEG  SCAN,  40  DB  TAPER,  20  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  62.  COMPARE  5  DEG  &  2.5  DEG  COMPENSATION 
0  DEG  SCAN,  40  DB  TAPE,  20  ELEMENTS 
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Compare  Compensation  Schemes 


BATED  AND  AMPLITUDE  ONLY  COMPENSATION 
DB  TAPER,  20  ELEMENTS 


Compare  Compensation  Schemes 
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Compare  Compensation  Schemes 


FIGURE  67.  COMPARE  UNCOMPENSATED  AND  AMPLITUDE  ONLY  COMPENSATION 
0  DEG  SCAN,  40  DB  TAPER,  40  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  68.  COMPARE  5  DEG  &  2.5  DEG  COMPENSATION 
0  DEG  SCAN,  40  DB  TAPER,  40  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  69.  COMPARE  COMPENSATED  AND  AMPLITUDE  ONLY  COMPENSATION 
45  DEG  SCAN,  40  DB  TAPER,  40  ELEMENTS 


Compare  Compensation  Schemes 


FIGURE  70.  COMPARE  5  DEG  &  2.5  DEG  COMPENSATION 
45  DEG  SCAN,  40  DB  TAPER,  40  ELEMEN] 


5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

Several  expected  relationships  were  evident  in  this  study.  Namely,  smaller 
r rays  dominated  by  edge  effects  and  are  therefore  "i  or:  ■-anti-/  rvpacted  by 

mutual  coupling.  Severe  amplitude  tapers  and  resulting  low  side  lobe  patterns 
are  difficult  to  achieve.  And  lastly,  scanning  compounds  these  effects.  From  a 
design  point  of  view,  the  following  guidelines  may  prove  useful. 

For  30  and  40  element  arrays,  mutual  coupling  may  be  ignored  for  side  lobe 
considerations.  Since  the  far-field  patterns  are  so  well  behaved,  the  inevita¬ 
ble  side  lobe  growth  at  large  scan  angles,  can  be  easily  dealt  with  by  specify¬ 
ing  side  lobe  performance  3  to  5  db  lower  than  necessary  at  broadside. 

The  20  element  array  just  crosses  the  arbitrary  3  db  growth  standard  at 
large  scan  angles,  and  severe  tapers.  The  compensation  techniques  provide  some 
measurable  improvement  up  to  30  db  Chebyshev  tapers,  however  near  40  db  (-35  to 
-38)  tapers  are  achievable  without  resorting  to  coupling  analyses  and 
compensation  techniques. 

The  10  element  arrays  are  dominated  by  coupling,  however  30  db  tapers  can 
be  achieved  using  amplitude  only  compensation  for  moderate  scan  volumes  <  30° . 
The  40  db  tapers  can  be  achieved  near  broadside,  but  are  not  achievable  for  45 
deg  scan  angles  using  any  of  the  compensation  schemes  in  this  report. 

It  is  interesting  to  note  that  the  amplitude  plus  5  degree  phase 
compensation  scheme  was  less  effective  than  the  amplitude  only  scheme.  The  2.5 
degree  scheme  was  approximately  the  same  as  the  amplitude  only  scheme,  therefore 
phase  compensation  techniques  provided  no  additional  side  lobe  performance  to 
justify  their  use  over  the  amplitude  only  technique  (for  the  array  geometries 
and  excitations  studied  in  the  report).  This  may  not  hold  true  for  greater  than 
40  db  Chebyshev  tapers. 
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Areas  of  further  study  could  either  extend  the  fortran  codes  in  this  study 
to  model  linear  arrays  of  strip  dipoles  with  specific  corporate  feed  structures, 
or  [rrp-r-r*  :f  ing  t  ?ore  complex  shaped  pattern  design 

problems.  Comparisons  of  amplitude  versus  phase  only  tapers,  as  well  as  the 
sensitivity  of  unequally  spaced  or  sparse  linear  arrays  to  coupling,  could 
easily  be  added  to  this  code.  If  planar  arrays  are  of  interest,  the 
King-Middleton  Five-Term  theory  would  be  required,  'which  would  entail 
considerable  modification  of  the  author's  programs.  Existing  method  of  moments 
codes  may  be  easier  to  apply,  in  many  instances,  than  the  King-Middleton  Theory, 
at  some  additional  computational  cost.  If  this  code  is  extended,  available 
symmetries  should  be  exploited  to  reduce  the  required  computer  resources. 
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APPENDIX  A1 


DERIVATION  OF  KING-MIDDLETON' S 
ZERO-ORDER,  TOO- TERM  THEORY  FOR 
DIPOLE  ANTENNA 


103 


Al.i  ISOLATED  DIPOLE 


Consider  the  antenna  geometry  shown  in  Figure  A.l.  The  following  condi¬ 
tions  are  imposed  upon  this  dipole: 


h  » 

a 

(A-la) 

(3. a  << 

1 

(A-lb) 

Iz{z)  = 

I.  (-z) 

(A-2a) 

A  (z)  » 

Az(-z) 

(A-2b) 

♦  (2)  = 

— <j>  ( — z ) 

(A-2c) 

q  (z)  = 

-q(-z) 

(A-2d) 

lz{±  h)  = 

0 

( A-3a) 

Z  s 
1 

0 

( A-3b) 

Equations  (A-la  &  b)  establishes,  in  a  general  sense,  what  is  sufficiently  thin. 

i 

Equations  (A-2a-d;  stare  one  symmetry  conditions  ./.at  ancle  to  .vlindrital 
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dipole  antennas.  Equation  ■; A-3a )  states  currents  vanish  at  the  end  of  the 
antenna  (which  can  be  achieved  in  practice  by  using  hemispherical  end  caps),  and 
■  A-3b>  ,ThicK  ctcr-'c  "cot  *\he  internal  impedance  of  a  perfect  conductor  ‘.s 
exactly  zero. 

Given  (A-3b),  the  vector  potential  must  satisfy  the  homogeneous  Helmholtz 
for  a  Z-aligned  dipole  (l-12a),  restated  here 


- +  0  * 

,  2  o 

I3z 


A  =  0 


( A-4 ) 


The  homoaeneous  solution  to  (A-4)  can  be  written  (15:48) 


p=a  ,  6  <  z  <  h 


A  = 
2 


f-j 

f 

\ 

— 

ci 

COS 

0O 

z  +  C, 

SIN0  z 

o 

V 

/ 

0 

-j 

/ 

Cx 

cos 

So 

2  -  C2 

SIN0  z 

0 

/ 

(A-5) 


p=»a  ,  -6  >  z  >  -h 


Substituting  (A-5)  into  the  Lorentz  condition  equation  (l-4c),  and  equating  the 
gap  voltage,  Vq  ,  to  the  scalar  potential  difference,  <j>(6)  -  4>(  — 5 ) ,  one  can  solve 
for  the  C2  coefficient.  Equation  (A. 5)  can  be  reformulated  for  a  "slice"  gener¬ 
ator  ( 6-0 ) . 

,  -h  <  z  <  h  .  (A-6 ) 

p=a 


-j 

A  (z)  =  —  C.  COS0  z  +  Vo  SIN0  z 

Z  i  O  -* —  O  * 

v  2 


From  retarded  potential  theory,  the  vector  potential  on  the  surface  of  a 
cylindrical  antenna  can  be  related  to  an  integral  over  all  currents  within  that 


cylinder  (14:76).  Further,  all  currents  within  that  cylinder  can  be  treated  as 
an  axial  current.  That  is 


A  (z) 


p-a 


z'  ) 


K ( z , z ' )  dz' 


where 

-30o  R 

K( z,z' )  =  e _ 

R 


R2  =  ( z-z'  )2  +  a2 


(A-7a) 


(A-7b) 

(A-7c ) 


K(z,z')  is  the  free  space  Green's  function,  and  Iz(z')  is  the  axial  distribution 
of  current,  and  z'  is  the  coordinate  of  source  points,  and  z  is  the  coordinate 
of  field  points.  By  equating  (A~7a)  with  the  general  solution  of  the  vector 
Helmholtz  equation  (A-6),  we  can  form  the  Hallen  integral  equation  (15:52-54). 


h 

Iz (z' )K(z,z' )dz' 

/ 

-h 


-j4n 


c,  cose  Z  +  —  SINS  z 

1  O  ^  O  1  1 

2 


,  -h  <  z  <  h 


(A-8) 


Since  no  closed  form  solution  to  (A-8)  exists.  King's  solution  to  the 
Hallen  integral  equation  is  based  on  a  careful  choice  of  current  distribution, 
and  expansion  functions  (15:51),  such  that 

h 

/ 

Iz  (z'  )K(z,z' )  dz'  =  j(z)g»(zi  .  (A-9) 

'  -h 
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Thus  reducing  the  integral  equation  co  an  algebraic  equation.  For  example,  if 
we  let 

J ( z )  a  iaF(z)  (A-lOa) 

be  an  assumed  distribution,  and  the  actual  distribution  can  be  expressed 

I(z')  =  I  H( z' )  .  (A-lOb) 

O 

We  could  express  I(z')  in  terms  of  J(z)  by  dividing  (A-lOb)  by  (A-lOa)  and 
solving  for  I(z)  or 

(I  H(  z '  n 

Hz' )  =  J(z)  — -  a  j(z)g(z,z')  .  (A-lOc) 

ll.F(z) 

The  function  g(z,z')  is  an  expansion  function  relating  actual  and  assumed  dis¬ 
tributions.  Define 

h 

f 

<|/(z)  &  g(z,z'  )K(  z,z'  )dz'  .  (A-lOd) 

/ 

-h 

Since  H(z')  is  not  known,  choosing  g(z,z')  exactly  is  not  likely,  therefore 
(A-lOd)  can  be  used  in  the  left  hand  side  of  (A-8). 

h  h 

/  f 

Hz' )K(z,z' )dz'  =  J(z)\j/(zj  +  ('Hz'  '-J(z)g(z,z'  ''jK(z,z'  Hz' 

J  ( A-ll ) 

-h  -h 
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The  difference  integral  cn  the  right  side  of  (A-ll)  reintegrates  the  lost  cur¬ 


rent  distribution  in  the  first  term.  If  q\z,z')  is  chosen  accurately  enough, 
the  finest  tarin  r»n  be  called  the  dominating  term,  and  the  difference  ;  ~ red 
will  be  a  negligable  correction.  The  equation  (A-ll)  can  be  solved  approxi¬ 
mately  for  the  expansion  function 


'lz(z) 


1 


h 


'  4n 

' 

ll 

Iz (z' )  K(z,z')  dz'  =  — - 

Vz) 

J  (  z  ) 

> 

p-a 


(A-12a) 


Due  to  the  R  in  the  denominator  of  the  integrand  (see  ( A— 7b ) ) ,  we  cam  claim  that 
the  vector  potential  (cn  the  surface)  is  governed  mainly  by  the  currents  near  or 
at  z'=z  (14:15-17).  Therefore  ^2(z)  is  mostly  real  and  constant  over  the 
antenna.  Therefore 

z )  =  4'(zr )  (A-12b) 

where 


\j/(zr)  is  the  ratio  of  vector  potential-to-axial  current  at  a  reference 
point,  zr 

Therefore  the  expansion  function  of  (A-12a)  becomes  an  expansion  coefficient. 


s 


1 


J(zr  ) 


h 

/ 

I  (z't  K  (z  ,z') 

7.  r 

-h 


dz ' 


(A-12c) 


The  reference  point  chosen  by  King  was  the  point  of  maximum  current  in  J(z). 


King  actually  starts  with  a  Hallen  difference  integral  equation  (which 
vanishes  at  z=»+h  by  definition)  (14:81-86). 


Az (z)-Az (h) 


3  0 


z-+h 


(A-13) 


Therefore 


'  -j4n 

V  ) 

o  •  1 

I  (z' )Kd (z,z' )dz'  - - 

ci 

cose0z  +  —  siNe0 1  z  |  +  u 

2  J 

j 

-h 


and 


(A-14a) 


U  *  -4n  Az (h)  » 


-n 

4n 


lz(z' )K(h,z' 5dz' 


(A-14b) 


-h 


K  (  z ,  z '  )  S  K( z,z' )  -  K(h,z') 


(A-14c) 


Applying  the  boundary  condition  of  (A-13),  the  Cx  coefficient  can  be  determined. 


(U  +  (V /2)  SlN0oh' 


cose  h 


(A-15) 


and  (A-lla)  can  be  re-written 


Iz (z' )Kd (z,z' )dz'  = 


-h 


-j4n 


n  cose,  h 


U(COS6,z-COSeoh)  +  —  SIN6o  ( h- 1  z  | ) 


( A-16 ) 


109 


To  choose  the  form  of  the  current  distribution,  break  the  left  side  cf  (A-14a) 
into  the  four  integrals. 


Iz (z' )Kd (z,z' )dz'  = 


-h 


rcose  r 


Iz(2') 


-h 


R  J 


dz'  +  j 


Iz(z') 


rsiNeoR) 


R  ) 


dz' 


-h 


V2' 


fcoseA) 


dz'  -  j 


-h 


It  (z' 


-h 


fsiNSA) 


dz' 


(A-17) 


where 


R^  =  (h-z')2  +  a2 

Reducing  each  integral  according  to  the  behavior  of  the  Kernal  (see  Figure  A.2), 
that  is 


cose  R 

O 

R 


=  ij/6(z-z'),  y  is  a  constant 


SIN0  R 

O 

*•» 


cos 


'6A 

>  2  „ 


(A-18a) 


(A-18b) 


110 


SIN(BNOTH)/BNOTH 


FIGURE  A.  2.  FUNCTIONAL  APPROXIMATION  OF 


Plugging  (A-18a  &  b)  into  ( A— 17 ) ,  we  can  reduce  to 


h 

( 

|  I  iz')K,(z(z')dz'  =  I  (zju/,  -  (COS  1/26  z  -  COS  1/26  h)u>,  (A-19a) 

z  a  z  o  0  2 

-h 


and 


% 


DP  2 


I  (z')COS 


l/260z‘ 


dz' 


(A-19b) 


Inserting  (A-19a  &  b)  into  (A-16)  we  can  write 


(COS  1/26  Z-COS  1/26  h)m, 

O  0  2 


-j4n 

*\c°S30h 


U(C0S6sh-C0S6oz) 


V 

+  ~  SIN6.  ( h- 1 z  | ) 
2 


(A-20) 


Solving  (A-17)  for  Iz(z),  we  obtain  King's  three-term  current  distribution 
(15:63-65)  of  the  form 


M*)  -  MMo.  +  Vo.  +  TDHoJ  (A"21a) 

where  ly ,  Tw  and  TQ  are  complex  coefficients  and  the  current  distribution 
includes  the  shifted  sinusoidal,  shifted  cosinusdal,  and  shifted  half-angle 
cosinusoidal  terms. 

Moz  A  SIN60  ( h- 1  z  |  )  (A-21b) 

A  ( COS 6  z  -  C0S6  h)  (A-21c) 

Hoz  a  (COS  1/2 60  z  -  COS  1/26, h  i  (A-21d) 
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Now,  inserting  this  three-term  form  in  (A-21a  through  d)  into  (A-16),  we 


I  (M  +  T  F  +  T  H  ) K. ( Z , Z ' ) dz '  = 

v  o  z  Uoz  D  o  z  d 


I(-H  f  , 

K  —  Iv(Me 
l  4rt 

-h 


+  T„F  +  T  H  )K(h,z')dzMF  +  — M 

Uoz  D  0  Z  r  <■ 


02  2 


(A-22a) 


where 


n  COS 6  h 

o  o 


Restating  the  form  of  the  integral  reduction  rules  in  (A-19a), 


(A-22b) 


Iz  (z'  )Kd  (z-z'  )dz'  =  Iz{z’)*1  -  Hoe^. 


2 


(A-23) 


Equation  (A-22a)  reduces  to 


^v^dR^oz  +  ^v^dl^ozj  +  ^v^U^dUR^oz  +  ^  v  ^d  U I  z  j 


+  T_T|>,n„H  +  jl 

1  v  D  dDR  o  z  J 


v  D  d  D 1 1 


rr-jn0iv  v 

ku~(™+Td*d)fo‘  +r”‘ 


(A-24 ) 


The  many  expansion  functions  are  defined  as  follows, 
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h 


\p 


J  Moi(z')K(h,z')dz' 
-h 


Foz(z')K(h,z')dz' 

-h 


h 


/ 

Hoz (z' )K(h,z' )dz' 
-h 


'P. 


dR 


-h 


KR{zr ,z' J-K^th.z' ) 


dz' 


Zr" 


)  ,  (3oh  <  n/2 

i-X/4 ,  0  h  >  n/2 


'P, 


d  I 


M  (z' 

O  Z  ' 


-h 


\[zz,z’)  -  Kj  ( h ,  z ' ) 


dz' 


z  -0 

r 


dUR 


1  'l 


*  (OJ 


F  (z' 

O  Z  ' 


Vzr*z')  -  Vh'Z'> 


dz' 


-h 


z  *0 

r 


r  i 


'p. 


dUX 


(2  )J 

os'  r  ' 


F  (z'  ) 

os'  ' 


K:  (z,  , z ' )  -  Kx  (h,z' 


dz' 


-h 


z  =0 

r 


(A-25a) 


(A-25b) 


(A-25c) 


(A-25d) 


(A-25e) 


(A-25f ) 


(A-25g) 
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(A-25h) 


h 


'  1  i 

/ 

H„(z') 

K(z  ,z' )  -  K(h,z' ) 

dz' 

P  <  7.  \ 

.  r  J 

V  l  V 

1 

Grouping  all  terms  in  (A-24)  by  current  distribution,  and  equating  the 
coefficients  to  zero,  we  form  three  equations  in  three  unknowns;  Iv ,  TQ  and  Ty . 


KV 


I  ill 
v  vdR 


=  0 


T  il/ 
u  vdun 


lCOSPohJ 


(U/  +  Til/  +  T  )  »  0 

U  D*D  ' 


TD^dDR  +  3^dl  +  ^U^dUI  +  ^D^dDI  ^  ^ 

Solving  (A-26a  through  c)  for  Iv ,  TQ  and  Tu , 


KV 


J2nVa 


2*dR  ^^RCos^oh 


(A-26a) 


(A-26b) 


(A-26c) 


(A-27a) 


'  ^dl^dUR0050^  - 
'*dD(*duRC0SBoh  ~  *u>  + 


(A-27b) 


*v*dD  -  H^dl 


^dD(^URC0SBoh  -  ♦«>  +  j%*dU 


I ' 


(A-27c) 


Having  determined  these  unknown  coefficients,  we  have  solved  for  the  cur¬ 
rents  on  a  thin,  gap  driven  dipole.  This  form,  as  expressed  in  <A-21a)  is 
King-Middle ton' i  zero-order,  three-term  theory.  To  obtain  even  greater  accuracy 


115 


in  the  formulation  of  dipole  current  distributions,  plug  the  zero-order  current 
expression  (A-21a)  into  ( A— 11 ) ,  for  both  I(z')  and  J(z)  in  the  difference  inte¬ 
gral  (15:63-65).  The  difference  ir.teor^i  can.  be  evaluated  to  formulate  a 
correction  current  distribution  which,  when  added  to  the  zero-order  solution 
creates  a  first  order  solution.  Likewise  a  second  order  solution  can  be 
formulated  by  plugging  the  first  order  solution  into  the  difference  integral 
term  of  (A-ll)  and  deriving  a  second  correction  current  distribution.  The 
second  order  current  distribution  is  just  the  sum  of  the  zero-order  solution  and 
both  correction  distributions.  King  and  Middleton  demonstrated  excellent 
agreement  between  their  second  order  theory,  and  the  precise  impedance  mea¬ 
surements  of  R.B.  Mack  (see  Figure  A-3).  The  second  order  theory  is  excessively 
cumbersome  when  applied  to  numerous,  coupled  elements,  therefore  it  has  not  been 
developed  with  any  rigor  here.  As  King  has  shown,  lower  order,  two-term 
theories,  when  properly  corrected  with  lumped  constants,  are  adequate  for  many 
array  environments  of  interest  (14:63-65). 

Al.2  TWO-TERM  THEORY 

For  the  isolated  dipole,  the  three-term  theory  is  particularly  needed  for 
antennas  near  antiresonance  (or  0oh=Nn).  However, this  report  will  focus  on 
resonant  antennas  (that  is  f30h-m/2)  for  which  King  has  demonstrated  that  a 
modified,  zero-order,  two-term  theory  is  sufficiently  accurate.  To  formulate 
this  simpler  theory  for  resonant  antennas,  first  note  that  current  form  in 
(A-21a)  is  indeterminate  at  resonance  ( 6ghn/2 ).  From  inspection  of  (A-27a),  we 
see  that  C0S6oh  approaches  zero  at  resonance,  and  so  lv  and  lz(z)  are  indeter¬ 
minate.  A  special  form  can  be  derived,  which  is  determinate  at  resonance,  by 
rearranging  (A-21a)  (15:63-65).  The  Moz  current  term  can  be  expanded  (via  the 
trigonometric  identity,  SIN(X+Y)  =  SINXCOSY  ±  C0SX3INY)  ,  and  the  COS6oh 
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factor  can  be  recombined  with  T..  and  T„  such  that 


— i  2nV 


I  (z)  = 


G  -  T.,  'F  -  TfH 

O  Z  J  C  Z  ~  0  z 

/ 


(A-28a) 


where 


Goz  (Z)  =  3IN6o  |  z  |  -  SIN@,h 


(A-28b) 


•  +  SIN0  h' 

U  o 

COS0  h 


(A-28C) 


=  T_ /COS 8  h 


(A-28d) 


The  equations  (A-28a  through  d)  are  determinate  at  resonance.  Notice  that  cur¬ 
rent  term  Moj,  has  been  replaced  by  Goz-  As  previously  stated,  a  two-term  formu¬ 
lation  is  sufficiently  accurate  for  far-field  and  driving  point  impedance  consi¬ 
derations.  To  achieve  this,  set 


Tu  +  Td 


in  (A-21a).  The  resulting  two-term  theory  can  be  expressed 


lz{z)  = 


-j2JlV. 


no*dP. 


G  +  T'  F 

1  O  Z  02 


(A-29a) 
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and 


SINS  L!SEC|3  h  +  . SIMS  h  - 

I  _  •>  U _ - _ _ 2 _ 1_ 

^  *u-*dUC0S6:h 


~.“u/ 


A1.3  TVJO-TERM  IMPEDANCE  FORMULAS 

To  formulate  the  admittance  of  the  isolated  dipole,  apply  Ohm's  Law  to  the 
gap,  and  evaluate  the  base  current  by  setting  Z=0  in  (A-29a  &  b).  That  is  the 
admittance  can  be  expressed 


Y  a 


VO) 

V 


-j2n  f 


T'f] 


z=0 


or 


Y 


-j2n 


'  \ 

-SIN0  h  +  T'U-COSB  h) 

o  o 

l 


(A-30) 


Figure  A. 3  compares  the  admittance  of  (A-30)  with  Mack's  measured  values.  The 
conductance,  G,  is  the  real  part  of  the  admittance,  Y,  and  the  subceptance,  B, 
is  the  imaginary  part. 
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Dipole  Admittance  Theories 
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FIGURE  A.  3.  COMPARISON  OF  KING-MI DDLETON  2ND  ORDER, 

&  ZERO-ORDER  THEORIES  WITH  MEASURED  DATA 


A1 . 4  TWO-TERM,  FAR- FI ELD  FORMULAS 

The  far-fields  of  the  isolated  dipole  can  be  formed  by  superposing  the  far- 
field  of  each  current  term.  3y  substituting  (A-39)  into  (1-19  and  b),  we  get 


E*  -  -  SIN© 

4rt 


r  -30op 

r 

'j2nv  r 

e 

l  r  J 

> 

A,*dR  ' 

j0  Z'COS© 


dz' 


-h 


or 


and 


K  = 


-V  f  -jeoR) 

-  e _ 

^dR  t  r 


f  (0,  6  h) 


( A-31a ) 
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-  1-CCS6  hCCS(f3  hCOS0)  COS0-SIN0  h  SIN(  (3  hCCSQ) 

l  5  ’  J  ’’ 

SIN©  COS© 

h 

60  SIN©  f  j(3  z '  COS© 

G  (0,6  h)  =  -  ( COS 6  z ' -COS6  h)e  °  dz' 

2 

-h 

SIN6. h  COS( 6. hCOS©)COS©  -  COS6.h  SIN( 6.  hCOS©) 

SIN©  COS© 

Ai . 5  DIPOLE  PAIR 

Consider  now  the  geometry  of  Figure  A.  4. 


( A-31c ) 


(A-31d) 
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The  same  symmetry  conditions  and  boundary  conditions  imposed  in  section  A.l  (A-l 
to  A-3)  are  imposed  on  both  dipoles  here.  If  we  form  the  Hallen  difference 
infect'd  !'A— for  coupled  diodes.  '•■;e  net  for  diode  J 

h 

(z' )K.  :d(z,z' )  +  Iz2  (z' )K:2d  (z,z' )jdz' 

/ 

-h 


=  K 


V,  SINg5(h-|z|  )  +  U.  (C0S6oz-C0S6oh)' 

.2 


and  for  dipole  2 


(A-32a) 


I2l(z')K21d!z,z')  f  I22(z')K22a(z,z') 


-h 


-  K 


V,  SINS  ( h-  I  z  1 )  +  U  (COS (3  _  -  COSS  h) 

c  O  1  1  i  0  2  O 

12 


where 


-jfloRn  -30oRilh 

Klld(z,z')  a  Klx (z,z' )-Kld (h,z' )  =  e  -  e 


Rir  «i 


lh 


K!2d 

(z,z' ) 

a 

K12(z,z' 

V 

K.  ,  .  '  z,z 

-  ‘  - 

-id 

K2!d 

( z, z' ) 

= 

Ki2d(z,z 

1  2 


^oR!2h 

3  -  e 


R1  2  R1  2  h 
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I,x (z' )K11 (h'2' )  +  I2J (Z' )K12 (h,z' ) 
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U, 
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Izl(Z')K21(h,Z')  +  Iz2(Z')K22(h'Z') 

J 

dz' 


-h 


(A-32e) 
( A-32f ) 


(A-32g) 


(A-32h) 


(A-32i) 


In  (A-32a  &  b),  we  have  two  simultaneous  integral  equations  to  solve  for 
the  current  distributions  on  each  dipole,  in  the  presence  of  the  other.  For  N 
coupled  elements,  (A-32a  &  b)  would  increase  to  N  simultaneous  integral  equa¬ 
tions,  and  the  lefthand  side  of  each  equation  would  integrate  over  N  current 
distributions,  Iti(z').  This  report  will  later  develop  a  matrix  formulation  for 
an  N  element  linear  array  of  parallel  dipoles.  At  this  point,  still  focusing  on 
the  two  coupled  antennas  of  this  section,  King  would  substitute  his  two  or  three 
-term  current  expansion  into  (A-32a-i),  depending  upon  the  actual  element  geome¬ 
try  (14:96-97,  133,  188,  207,  237-238,  272-280). 

Let  us  restate  the  two-term  current,  for  the  ith  dipole  as  follows 

I2l (z)  S  Iui(z)  +  I.JX(Z)  =  iA.n_  -  3  r  1  A— 33 ) 
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and  M  (z)  &  F  (z)  are  defined  in  (A-21b)  and  (a-21c)  respectively.  Substitut- 

O  Z  0  2 

ing  (A-33)  into  (A-32a-i),  numerous  integral  terms  over  Iui(z)  and  Iyi(z)  must 
be  reduced  to  algebraic  equation^ ,  'be  in tec re 3  ’-eduction  rules  vary  now  with 
the  spacing  ,b,  between  the  coupled  elements.  In  general,  for  30b  >  1 


)dZ' 


-  J  I  ,  (  Z ) 

yKidu  uk 


(A-34a) 


and  for  0  b  <  1,  the  reduction  rules  are 

o  r 


eob  >  i 


W2'^^2'2'^2' 


f  jA. 

- —  U/  I  (  Z  ) 
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l  B,  ) 


(A-34b) 


Iui(zMKKid(Z,2')d2'  =  —  *Kldu  Iuk(*) 


(A-34c) 


,i(2')KKidR(2'2')d2'  =  ~  *KidR  1,y{Z) 


6b  <  1 


(A-34d) 


^i^'^ldl^'2'^2'  =  -  *KidI  W2) 

l  B.J 


( A-34e ) 
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Additional  reduction  rules  needed  to  evaluate  U,  and  U_  in  (A-32h  &  i)  are, 


I  i,  z '  )K  .  (h,z'  )dz'  =  3  v, 

u  1  xi  i  k 


{ A — 34  f ) 


^1(2')^.  (h,z')dz'  =  jA.^ixj 


(A-34g) 


Commonly,  3Qb  <  1  is  only  true  for  i=k,  therefore  3ofc£0os  <<  1.  This  is  not 
generally  true,  but  element  spacings  will  be  chosen  to  meet  this  condition 
through  out  this  report. 

For  dipole  1,  (A-32a)  becomes 


Iul(z')  +  Iwl<z')  Klld(z,z')  +  Iu2(z')+Iv2(z')  K12d(z,z')  dz' 


M  +  KU,  F 

O  Z  1  O  2 


(A-35a) 


After  reducing  integral  terms,  we  have 


^lldu  ^ul^2*  +  ~  ^lldR  ^vl^2^  +  ^lldl  ^ul^2^ 

J3^  /  'A^ )  , 


I  .  I  J'-  2  i  u  -  -  _  I 
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(A-35b) 


M  +  KU,  F 

oz  i  ~  z 


Grouping  terms  for  dipole  1,  we  can  write 


-  \l/  -f” 

ylldu 


rjM  (Z2)  fjA,] 

-  *lldl  +  ~  *12du  +  - 

>  )  UB1  J  l  B,  J 


*»dv  -  rai  Fo 


(A-36a) 


d,  I  (z)  =  -  M 

^iidK  ^  oz 


(A-36b) 


Likewise  for  dipole  2,  equations  like  (A-36a  and  b)  can  be  written 


^2  IdR  + 


IT”-  *  ]*»„  '.i'*'  -  F°z 


(A-36c) 


U*2«  Iv2(2)  =  2  "" 


(A-36d) 


In  (A-36a  through  d),  we  have  four  equations  in  four  unknowns,  A1 ,  A2 ,  Bx ,  and 
B  .  We  can  then  write  for  the  Kth  dipole 


I  „(z) 

v  k 


(A-37a) 


B0  fjMf 

I’^idu  x  —  ^,'H 

V  IbJI 
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-  K 


(A-37b) 


r-jn,  2 
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— -  z 

+  Mklu 

,  4n  i=l 

where 


6 


ik 


a. 


0 ,  for  i^K 
.1 ,  for  i=K 


By  cancelling  out  the  current  distribution  on  either  side  of  (A-37b),  and 
grouping  terms  by  coefficients,  A.  &  B.  ,  we  can  write 


2 

l 

i-1 


B. 


*kiduC°S<* 


h  -  'P, 


k  i  u 


2 

• 

f  \ 

-  j  2 

*kiv- 

^kidM  ^  1— &i  k  )  +  1  ^kidl  ^ik 

COS  3  h 

O 

i-1 

. 

. 

(A-37C) 


From  (A-37a)  we  cam  note  that  the  A^s,  and  therefore  the  shifted  sinusoi¬ 
dal  distribution  term,  is  determined  entirely  from  the  gap  voltage,  and  the 
dipole  dimensions  (embodied  in  'l'kkdR).  In  (A-37c),  we  have  a  set  of  two  linear 
algebraic  equations  that  can  be  solved  for  two  unknowns,  Bx  and  B2 ,  in  terms  of 
Aj  and  Aj .  Upon  solving  these  linear  equations,  we  have  determined  the  currents 
on  both  dipoles,  and  can  compute  self,  mutual  and  driving  point  impedances,  in 
the  manner  presented  in  Section  one  (1-24  and  1-25).  The  far-fields  from  these 
coupled  dipoles  cam  be  determined  in  a  similar  manner  to  the  isolated  dipole. 

For  coupled  dipoles,  one  must  superpose  the  fields  of  each  current  distribution 
(see  equations  A-31a  through  d)  from  each  element  (see  equations  1-26  and  1-27) 
given  arbitrary,  complex  voltage  ’.-.’eights  'par*:  f  We  ~n  ccmcute  the  far- 

field 
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h 

j<4u,  2  \(  —  j  (3  (  r-z'COS©) 

Eg  =  -  SIN©  Z  I..(z')  e _  dz'  .  (A-38a) 

4n  i=l  “  r 


Substituting  (A-33),  and  separating  the  amplitude  and  phase  dependent  portions 
of  the  Kernal,  we  have 


where  F  (0,6h)  is  defined  by  ( 1— 28b)  and  G  (0,8  h)  is  defined  by  (A-31d). 
Notice  that  (A-38c)  cannot  be  separated  into  an  array  factor  and  a  unique  field 
factor  as  in  (l-28a). 
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APPENDIX  A2 

EQUATIONS  FOR  KING-MIDDLETON ' S 
ZERO-ORDER,  TWO- TERN  THEORY  FOR 
LINEAR  ARRAYS  OF  DIPOLES 


128 


Equation  (A-37c)  could  easily  be  extended  to  a  set  of  N  linear  algebraic 
equations,  for  N  coupled  dipoles.  As  N  increases,  matrix  notation  becomes  more 
efficient.  Assuming  ’■he  two- term.  current  form  of  (A-33),  we  can  express  the 
gap  driving  voltages  (which  are  complex)  in  vector  form. 


V  * 


V. 


u,  J 


(A-39) 


The  A.  coefficients  can  be  determined  from  the  vector  multiplication 


A  s 


r  a. 


1 


-jK'V 


(A-40a) 


K'  - 


where  K  is  defined  by  (A-22b) .  The  vector  of  B.  coefficients 


B  ■ 


(A-40b) 


is  related  to  the  A.  vector  by  a  set  of  linear  equations  like  (A-37c).  For  N 
elements,  these  equations  can  be  represented  in  matrix  form 


(A-41 ) 


or  by  matrix  inversion 
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4>  A 


(A-42) 


where 


[M 


.*'•]■  r +iiu. 


♦ki.  *  *kidu  COSBoh  - 


(A-43a) 


♦kiv  a  ^Kiv  -  ^KiDvcose0H(i-5IK)  -  J  *kidicosp0h  sIK 


(A-43b) 


The  admittance  matrix  of  self  and  mutual  admittances  of  the  N  elements  can 
also  be  expressed  in  matrix  notation.  Starting  with  a  vector  of  current  distri¬ 
butions,  that  is 


Vz)  "  f  xi(z)  j  “  [  *u  j"1  [  K  pAF0I!  +  jAM0B 


(A-44a) 


IN(z) 


By  substituting  (A-40a)  into  (A-44a),  we  get 


I  ( z )  ■  $ 

8  '  TU 


<j>  K'F  +  K'  I  M  V 

^  02  O  Z 


(A~44b) 


To  evaluate  driving  point  properties,  evaluate  (A-44b)  for  z=0,  or 
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(A-44c) 


where 

I  is  a  vector  of  base  currents 

O 

[I]  is  a  NXN  identity  matrix. 

The  admittance  matrix  is  then  derived  by  the  matrix  expression  inside  the  outer 
brackets. 
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(A-45) 


The  impedance  matrix  is  just  the  inverse  of  the  admittance  matrix. 


(A-46) 


It  is  interesting  to  notice  that  the  admittance  matrix  is  solely  a  function 
of  array  geometry  and  the  current  distributions.  The  driving  point  characteris¬ 
tics  of  an  antenna  are  determined  by  the  driving  voltage  (and  resultant  cur¬ 
rents)  of  all  antennas. 
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APPENDIX  B 

UNCOMPENSATED  MODEL  RESULTS:  10  ELEMENT  ARRAY 

20  ELEMENT  ARRAY 
30  ELEMENT  ARRAY 
40  ELEMENT  ARRAY 
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APPENDIX  C 

COMPENSATION  RESULTS:  10  ELEMENT  ARRAY 

20  ELEMENT  ARRAY 
40  ELEMENT  ARRAY 
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